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The  problem  of  vehicle  mobility  is  a complex  one  in  that  a rigorous  analyti- 
cal treatment  of  the  mechanics  Involved  between  vehicle  characteristics, 
such  as  geometry,  size,  and  driving  forces,  smd  the  properties  and  responses 
of  soil  that  supports  a moving  vehicle  has  not  been  developed.  In  general. 
Immobilization  problems  occur  as  a result  of  the  loss  of  or  excessive  demand 
for  traction  created  by  a combination  of  slnkage,  slope,  forces  applied,  ob- 
stacles, and  environmental  conditions.  Immobilization  problems  associated 
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20.  ABSTRACT  (Continued) 

with  poor  traction  due  to  weak  coils  or  slippery  surfaces  are  not  limited  to 
rigid  or  pneumatic -tired  wheels,  but  also  affect  powered  track  wheels  with 
grousers.  Two  approaches  have  been  followed  in  attempting  to  solve  the  soil- 
wheel  interaction  problem:  an  nnnlyt.ical  approach  based  on  the  elastic  solu- 

tion of  a plane  strain  problem,  and  an  experimental  approach  based  on  the 
relationship  between  the  shear  and  normal  stresses  that  may  occur  within  the 
vicinity  of  or  at  the  interface  of  the  soil  and  the  wheel  of  a moving  vehicle. 
The  analytical  solution  is  based  on  the  assumption  that  the  stresses  within 
the  soil  are  the  results  of  the  tangential  and  radial  stresses  created  by  a 
wheel  partially  embedded  in  soil  with  an  Imoendlng  rotation.  The  Airy  stress 
function  was  used  in  representing  the  stres,;es  within  the  soil  In  terms  of 
analytic  functions.  The  Schwarz-Christoffe  . equation  was  used  to  transform 
the  geometry  and  the  boundary  condition  of  the  region  beneath  the  wheel  and 
to  match  them  with  the  stress  functions.  The  Cauchy  Integral  equation  was 
applied  on  the  transformed  boundary  conditions  to  obtain  the  shear  and  normal 
stress  at  any  point  within  the  region  of  the  soil-wheel  system.  A computer 
program  for  reducing  the  results  and  obtaining  numerical  values  of  the 
stresses  at  any  point  within  the  vicinity  of  the  wheel  was  also  written.  It 
is  believed  that  the  analytical  solution  developed  will  permit  the  evaluation 
of  stresses  within  the  soil  beneath  a wheel  that  result  from  various  combina- 
tions of  radial  and  tangential  stresses.  The  experimental  investigation  was 
designed  to  Investigate  the  shear  stresses  and  traction  forces  that  may  exist 
between  a model  rigid  wheel  or  tire  wheel  and  the  supporting  soil.  CH  mate- 
rial (Vicksburg  buckshot  clay)  was  compacted  and  tested  In  an  annular  torsion 
shear  machine.  Four  types  of  specimens  were  tested  in  this  study:  the  first 

series  of  tests  consisted  of  shearing  homogeneous  soil  specimens;  in  the 
second  series,  the  soil  was  sheared  against  smooth  rubber;  in  the  third 
aeries,  the  soil  was  sheared  against  polished  stainless  steel;  anu  in  the 
fourth  series,  a nonhomogeneous  soil  specimen,  with  upper  and  lower  halves 
having  different  water  contents,  was  sheared.  The  initial  normal  stresses 
used  in  the  teat  program  were  5,  15,  and  30  psl,  and  all  specimens  were 
sheared  under  constant-volume  conditions  at  rates  of  shear  deformations  of 
0.002  , 0.2,  and  2.0  in. /min.  Test  results  showed  that  nonhomogeneous  soil 
specimens  with  upper  and  lower  halves  having  water  contents  of  26  and 
16  percent,  respectively,  were  stronger  than  homogenous  soil  specimens  with  a 
uniform  water  content  of  26  percent  and  that  the  strengths  of  both  types  of 
soil  specimens  were  higher  than  those  at  the  Interface  of  smooth  rubber  or 
polished  stainless  steel  and  soil.  Test  results  showed  that  the  strength  of 
nonhomogeneous  soil  specimens  was  higher  than  that  of  homogeneous  soil  speci- 
mens with  uniform  water  content,  and  the  strengths  of  both  types  of  soil 
specimens  were  higher  than  those  at  the  Interface  of  smooth  rubber  or  pol- 
ished stainless  steel  on  soil.  The  results  also  showed  that  the  shear  resis- 
tance developed  between  soil  and  rubber  is  much  higher  than  that  developed 
between  soil  and  polished  steel.  The  shear  stress  at  failure  for  all  test 
series  showed  an  increase  with  increasing  initial  normal  stress  and  increas- 
ing rate  of  shear  deformation.  The  shear  stresses  at  the  interfaces  of 
smooth  rubber  and  soil  and  also  of  polished  steel  and  soil  dropped  sharply 
when  the  plane  of  contact  was  covered  with  a film  of  water.  The  peak  shear 
stress  obtained  from  each  test  series  was  plotted  as  ti  function  of  the  cor- 
responding normal  stress  in  order  to  facilitate  its  use  in  analytical 
solutions. 
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PREFACE 


The  study  reported  herein  was  conducted  at  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES)  under  the  sponsorship  of  the  Office, 
Chief  of  Engineers,  Directorate  of  Military  Engineering,  Project  No. 
UA061102B52E,  "Research  in  Military  Engineering  and  Construction," 

Task  01,  "Environmental  Characterization  System  Research,"  Work  Unit  013, 
"Field  Determination  of  Shear  Strength  Parameters  and  Traction  in  Near- 
Surface  Soils,"  Work  Unit  Accession  No.  DA  of  8156. 

The  analytical  and  experimental  studies  described  in  this  report 
were  performed  during  the  period  September  1970  to  June  1971  by  Dr.  M.  M. 
Al-Hussaini,  Project  Engineer,  and  Mr,  P.  A.  Gilbert,  both  of  the  Em- 
bankment and  Foundation  Branch,  Soils  and  Pavements  Laboratory  (SiPL). 
This  report  was  prepared  by  Dr.  Al-Hussaini,  assisted  by  Mr.  Gilbert, 
under  the  immediate  supervision  of  Mr.  J.  R.  Compton,  Chief,  Embankment 
and  Foundation  Branch,  and  the  general  supervision  of  Mr.  J.  P.  Sale, 
Chief,  and  Mr.  S.  J.  Johnson,  Special  Assistant,  S&PL.  Useful  sugges- 
tions and  comments  by  Dr.  G.  Y.  Baladi,  S&PL,  are  appreciated. 

Directors  of  WES  during  the  course  of  the  investigation  and  prepa- 
ration and  publication  of  this  report  were  BG  E.  D.  Peixotto,  CE,  and 
COL  G.  H.  Hilt,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Sl) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  In  this  report  can  be  con- 
verted to  metric  (Sl)  units  as  follows: 


Multiply 


To  Obtain 


inches 

square  Inches 

pounds  (mass) 

pounds  (force) 

poimds  (force)  per 
square  inch 

pounds  (mass)  per 
cubic  foot 

Inches  per  minute 

Fahrenheit  degrees 


O.O25U 
0.00061*516 
0.1*535921* 
1*.  1*1*8222 


meters 

square  meters 

kilograms 

newtons 


6891*.  757  pascals 


16.0181*6 

0.0251* 

5/9 


kilograms  per  cubic  meter 
meters  per  minute 
Celsius  degrees  or  Kelvins* 


• To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read- 
ings, use  the  lollowing  formula:  C * (5/9)(F  - 32)  . To  obtain  Kel- 
vin (K)  readings,  use:  K ■ (5/9)(F  - 32)  + 273.15  • 
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STRESSES  AND  SHEARING  RESISTANCE  IN 


SOIL  BENEATH  A RIGID  WHEEL 


PART  I : INTRODUCTION 

Background 

1.  The  prediction  of  the  stress  distribution  beneath  a wheel  is 
important  in  soil-wheel  interaction  studies  since  the  magnitudes  of 
radial  and  tangential  forces  generated  by  the  vehicle  are  the  governing 
factors  in  the  stability  of  the  supporting  soil,  which  the  ultimate  mo- 
bilization of  a vehicle  depends  upon.  Two  approaches  are  known  to  be 
used  in  prediction  of  stresses  beneath  wheels  or  tracks  of  military 
vehicles: 

a.  The  analytical  approach  in  which  the  soil  is  Idealized  and 
assumed  to  have  elastic^  or  plastic  properties. ^ 

b.  The  experimental  approach  in  which  pressure  or  stress 
cells  are  used  to  measure  the  vertical  component  of  stress 
at  the  contact  area  of  the  wheel3,l^  or  beneath  the  soil 
surface . 1 

Analytical  approach 

2.  The  analytical  analysis  of  stress  distribution  within  the  soil 
beneath  a military  vehicle  is  mainly  based  on  Boussinesq's  theory,^ 
which  is  derived  to  evaluate  the  stress  distribution  beneath  a concen- 
trated load  acting  vertically  on  the  horizontal  plane  of  a semi-infinite, 
elastic,  homogeneous  material.  According  to  Boussinesq's  theory,  the 
vertical  stress  a at  any  point  within  the  medium  can  be  expressed  as: 


where 

P * Intensity  of  the  applied  load 
r ■ radial  distance  from  the  load 

y ■ depth  of  the  point  below  the  surface  as  shown  in  fig.  la 
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Equation  1 has  been  Integrated  by  Egorov"  to  obtain  the  stress  and  dis- 
placement beneath  a uniformly  loaded  circular  area.  The  special  case 


of  vertical  stress  a 

y 

be  expressed  as  follows: 


beneath  the  center  line  of  the  loaded  area  can 


Cy-Pn. 


W * J 


a > radius  of  the  loaded  area  as  shown  in  fig.  lb 
3.  The  moat  attractive  aspect  of  equations  1 and  2 is  that  the 
vertical  stress  at  any  point  is  only  a function  of  the  shape  and  the 
Intensity  of  the  applied  load  and  also  the  distance  of  the  point  from 
the  surface  but  not  of  the  elastic  properties  of  the  medium.  However, 
the  Bousslnesq  and  Egorov  equations  were  derived  for  specific  conditions, 
and  unless  their  basic  assumptions  are  fully  satisfied,  their  random 
application  in  mobility  problems  may  lead  to  erroneous  results.  Perti- 
nent points  that  must  be  kept  in  mind  before  applying  equation  1 or  2 in 
evaluation  of  stress  distribution  beneath  a wheel  are; 

a.  The  projection  of  the  contact  area  between  the  wheel  ;nd 
the  soil  is  neither  a point  nor  a circle  but  rather  an 
ellipse. 

b.  The  theories  were  derived  for  a semi-infinite  medium  with 
a horizontal  plane  surface,  and  they  do  not  take  into  con- 
sideration the  geometry  at  the  soil  surface  created  by  the 
moving  wheel. 

c.  The  theories  do  not  account  for  the  shear  stresses  on  the 
contact  area  between  the  wheel  and  soil,  which  are  re- 
sponsible for  the  traction  resistance  of  a moving  vehicle. 

d.  The  applied  stresses  are  usually  in  the  radial  and  tangen- 
tial direction  of  the  wheel-soil  contact  area  and  are  not 
perpendicular  to  the  horizontal  surface  of  soil. 

Experimental  approach 

U.  Because  of  the  limited  application  of  the  Bousslnesq  and 
Egorov  theories  in  mobility  study,  more  emphasis  was  placed  on  experi- 
mental stress  analysis.  Indeed,  many  experimental  devices,  such  as  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  pressure  cell. 


have  been  used  in  measuring  vertical  stresses,  at  the  surface  and 
within  the  soli  mass,  induced  by  a stationary  or  moving  vehicle.  Pres- 
sure cells  used  are  either  operated  by  fluid,  such  as  mercury,  or  are 
electronic  strain  gages;  they  vary  in  shape,  size,  and  accuracy,  and 
they  differ  even  in  the  way  they  are  placed  In  soil  to  obtain  better 

I4 

resolution.  The  description  and  the  advantages  of  each  type  of  pres- 
sure cell  are  beyond  the  scope  of  this  report.  However,  regardless  of 
the  type  pressure  cell  used,  there  Is  always  physical  limitation  in- 
herent In  this  branch  of  study  that  makes  pressure  cell  readings  highly 

0 

inconsistent.  According  to  Knight,  inconsistency  in  the  measured  data 
can  be  attributed  to  any  one  or  a combination  of  several  factors,  such 
as  variable  path  of  vehicle  due  to  guide-channel  movement  in  the  soft 
soil,  pressure  cell  movement  during  the  response  period,  stress  concen- 
tration or  dispersion  due  to  variable  soil  conditions  above  and  adjacent 
to  the  cell,  or  improper  seating  of  the  pressure  cell.  Other  such  fac- 
tors are  the  size  of  the  cell  compared  to  the  maximum  size  of  soil  par- 
ticles and  the  sensitivity  of  pressure  cell  reading  to  lateral  and  shear 
stresses  that  cannot  be  measured  by  the  existing  pressure  cells.  Thus 
it  appears  that  In  addition  to  the  enormous  cost  and  time  involved  in 
experimental  stress  analysis  using  pressure  cells,  it  is  almost  impos- 
sible to  obtain  a precise  measurement  with  the  presently  available 
instrumentation. 

Purpose 

5.  From  the  above  discussion,  it  can  be  stated  that  current  ana- 
lytical and  experimental  techniques  are  inadequate  to  provide  data  input 
for  the  design  of  Improved  military  vehicles  from  the  standpoint  of  their 
ability  to  travel  across  natural  terrain.  Such  an  objective  cannot  be 
achieved  unless  the  stress  distribution  at  the  surface  and  within  the 
soil  as  it  is  deformed  by  the  moving  wheel  is  thoroughly  understood. 

This  objective  represents  in  part  the  purpose  of  this  report. 

6.  The  primary  purposes  of  this  investigation  were: 
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a.  To  present  a new  theoretical  approach  for  determining  the 
magnitudes  and  the  orientations  of  normal  and  iheur 
stresses  beneath  a circular  depression  when  acted  upon  by 
uniform  radial  and  tangential  stresses. 

b.  To  generalize  the  solution  to  Include  the  evaluation  of 
normal  and  shear  stresses  under  any  combination  of  radial 
and  tsingentlal  stresses. 

c. .  To  use  the  analytical  solution  in  solving  practical  prob- 

lems of  mobility  studies,  such  as  towed  and  power  wheels, 
and  to  provide  theoretical  reasons  to  explain  their  ob- 
served performance  variation. 

d.  To  stimulate  interest  in  the  development  of  a more  closed- 
form  solution  based  on  broadly  accepted  stress-strain  re- 
lations that  can  be  used  more  systematically  and  logically 
in  terrain-vehicle  mechanics  as  well  as  applied  soil 
mechemlcs. 

To  investigate  the  physical  properties  of  typical  soil  and 
to  estimate  the  relationships  between  the  normal  suid  shear 
stresses  that  can  be  expected  at  the  interface  between  the 
soil  and  a wheel  of  a moving  vehicle. 

f_.  To  study  the  influence  of  different  factors,  such  as  water 
content,  rate  of  shear  deformation,  density,  sxirface  wet- 
ness, and  initial  stress  levels,  on  the  relationship  be- 
tween the  shear  and  normal  stresses  that  occur  with  the 
soil-wheel  system. 


Scope 


7.  Theoretical  solutions  were  developed  in  this  study  to  investi- 
gate the  distribution  of  normal  and  shear  stresses  within  a semi- 
infinite  mass,  the  upper  boundary  of  which  consisted  of  a segment  of 
circular  arc  connected  with  two  straight  lines  extended  to  Infinity,  in 
an  effort  to  generate  a reasonable  analytical  solution  to  the  wheel-soil 
intercu:tion  problem.  The  circular  arc,  which  simulates  the  area  of  con- 
tact between  the  moving  wheel  (i.e.,  rigid  wheel  or  pneumatic  tire  with 
inflation  pressure  close  to  average  contact  pressure),  is  acted  upon  by 
a uniform  radial  or  tangential  or  a combination  of  both  stresses.  Re- 
sults of  normal  and  shear  stresses  within  the  framework  of  the  proposed 
wheel-soil  system  are  presented  as  influence  tables;  thus,  their  usage 
will  not  be  restricted  to  specific  cases.  Available  analytical  and  ex- 
perimental data  provided  by  WES  €Uid  other  investigators  were  also 
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correlated  with  the  new  solutions  whenever  possible. 

8.  The  experimental  study  was  limited  to  CH  material,  compacted 
at  water  contents  of  26  and  30  percent,  which  correspond  to  dry  densi- 
ties of  95  and  92  pcf,*  respectively.  Because  laboratory  tests  were 
Intended  to  simulate  soil-wheel  field  conditions  as  closely  as  possible. 
It  was  decided  to  shear  the  soil  specimens  under  constant  volume  for 
large  jhear  deformations  and  at  rates  of  strain  higher  than  those  com- 
monly used  in  soil  testing.  The  WES  19^7  annular  shear  apparatus  was 
chosen  for  the  test  program  after  equipping  it  with  a constant-volume 
loading  system  and  automatic  recording  system.  The  experimental  study 
was  limited  to  shearing  of  (a)  homogeneous  soil  specimens,  (b)  soil 
against  tire  rubber,  (c)  soil  against  steel,  and  (d)  soil  specimens 
Inhomogeneous  with  respect  to  water.  The  initial  normal  stresses  im- 
posed on  soil  specimens  at  the  beginning  of  shear  were  limited  to  5,  15, 
and  30  psl,  which  cure  comparable  to  those  normally  encountered  under 
moving  vehicles.  The  rates  of  shear  deformation  adopted  In  the  test 
program  were  0.002,  0.200,  and  2.00  in. /min. 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Sl)  units  Is  presented  on  page  5. 
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PART  II:  STRESS  DISTRIBUTION  BENEATH  A UNIFORMLY 

LOADED  CIRCULAR  SEGMENT  IN  A SEMI-INFINITE  MASS 


9>  WES  was  assigned  the  task  of  conducting  experimental  tests 
simulating  the  movement  of  wheels  on  soft  soil  in  an  attempt  to  find 
some  meaningful  way  to  Improve  design  of  military  vehicles  for  mobility 
in  rugged  terrain.  This  task  was  an  extension  of  trafflcablllty  of  soil 
studies  that  have  been  taking  place  at  WES  for  many  years.  As  a result 
of  the  experimental  study,  an  effort  was  made  to  explore  and  determine 
the  Interdependent  stress  relationship  that  exists  between  the  wheel  and 
the  supporting  soil. 


Statement  of  the  Problem 


10.  In  this  study,  attention  was  focused  toward  the  determination 
of  stress  within  the  soil  mass  from  the  knowledge  of  contact  stresses 
beneath  the  moving  wheel.  The  first  step  consisted  of  a study  of  uni- 
form radial  and  tangential  contact  stresses  that  will  be  expanded  later 
on  to  Include  any  form  or  combination  of  contact  stresses  desired.  Al- 
lowance will  be  made  to  extend  the  solution  for  the  determination  of  a 
strain  field  within  the  soil-wheel  system  in  order  to  provide  better 
understanding  of  the  soil  deformation  pattern  beneath  the  wheel,  such  as 
slnkage,  bulldozing,  etc.,  on  a more  rational  and  more  exact  basis. 
Assumptions 

11.  As  an  introduction  to  the  analytical  solution,  the  following 
assumptions  were  made: 

a.  The  contact  surface  between  the  soil  and  the  wheel  has  a 
circular  curvature. 

b.  The  soil  beneath  the  wheel  is  a modified  seml-lnflnlte  mass. 

£.  The  soil  is  weightless  and  isotropic. 

d.  The  contact  stresses  between  the  wheel  and  soil  are  uni- 
form radial  and  tangential  stresses. 

e.  A state  of  plane  strain  conditions  exists  within  the 
framework  of  the  wheel-soil  system. 

12.  The  assumptions  stated  above  are  Justified  as  follows. 
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The  deformation  of  a rigid  wheel  or  pneumatic  tire  with  medium  to  high 
inflation  pressure  is  not  appreciable  compared  to  the  deformation  of 
soft  soil.  As  a result,  the  deformed  shape  of  soil  tends  to  take  the 
shape  of  the  wheel,  which  Justifies  the  first  assumption.  The  second 
assumption  is  Justified  because  stresses  are  significant  only  within  the 
vicinity  of  the  wheel  and  diminish  very  rapidly  with  increasing  distance 
from  the  surface.  The  soil  region  is  called  a modified  seml-lnflnlte 
mass  because  of  the  presence  of  circular  depression  (see  fig.  2)  in  the 
half-space.  The  assumption  of  weightless  material  is  Justified  simply 
because  stress  due  to  wheel  load  is  concentrated  at  regions  close  to  the 
surface,  whereas  effect  of  stress  due  to  the  weight  of  soil  is  immate- 
rial. However,  the  effect  of  the  weight  of  material  becomes  more  slgnl- 
cant  at  great  depth. 

13.  The  assumption  that  the  contact  radial  stress  is  assumed  to 

3 I* 

be  uniform  is  based  on  actual  WES  test  data  * that  show  that  the  ver- 
tical stresses  along  the  central  portion  of  the  contact  area  in  the  di- 
rection of  motion  are  fairly  constant  and  proportional  to  the  inflation 
pressure.  Since  the  normal  contact  pressure  is  constant  and  the  coef- 
ficient of  friction  between  a given  soil  and  wheel  is  constant,  then  the 
resulting  sheai  stress  has  to  be  uniform.  Thus,  the  assumption  of  uni- 
form tangential  stress  is  also  Justified. 

it).  In  defense  of  the  last  assumption,  it  could  be  said  that  ver- 
tical and  horizontal  components  are  at  their  siaximum  in  the  direction  of 
the  applied  stress,  i.e.,  in  a plane  parallel  to  the  vehicle  motion. 
However,  the  deformation  in  planes  transverse  to  vehicle  motion,  even  at 
the  surface,  is  almost  negligible  compared  to  that  in  planes  in  the  di- 
rection of  motion.  Thus,  the  existence  of  plane  strain  conditions 
within  the  soil-wheel  system  is  not  far  from  reality. 

Theoretical  considerations 

15.  In  the  previous  section,  it  was  shown  that  the  wheel-soil  in- 
teraction problem  can  be  closely  approximated  as  plane  strain  cases. 
Thus,  equations  for  plane  strain  theory  of  elasticity  are  applicable. 
However,  before  an  analytical  solution  to  plane  strain  problems  can  be 
obtained,  the  conditions  of  equilibrium,  boundary,  and  compatibility 
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must  be  satisfied.  These  three  conditions  constitute  the  necessary  and 
sufficient  conditions  to  satisfy  the  solution  of  an  elasticity  problem. 

l6.  Equilibrium  conditions.  Consider  a soil  element,  at  point  A 
beneath  a wheel  (see  fig.  2),  to  be  under  plane  strain  conditions.  This 
element  (see  fig.  3a)  is  considered  to  be  so  small  that  loading  is  es- 
sentially uniform  and  representable  in  the  first  approximation  by  a 
single  load  in  the  center  of  the  plane  of  application.  Since  the  mate- 
rial is  considered  to  be  weightless,  the  equilibrium  condition  can  be 
obtained  by  summing  forces  in  the  x and  y directions,  resulting  in 
the  following: 

3o  3t 
3o  3t 


a and 

X 


normal  stresses  in  the  x and  y directions, 
respectively 

shear  stress 


17.  Boundary  conditions.  The  stress  components  vary  throughout 
the  soil  mass,  and  equilibrium  conditions  must  be  satisfied  at  any  point 
within  the  soil  mass,  including  the  boundary.  In  this  respect,  the 
boundary  stresses  are  considered  to  be  a continuation  of  the  internal 
stress  distribution.  Assuming  BC  is  any  small  element  that  coincides 
with  the  boundary,  then  the  boundary  conditions  from  fig.  3b  are: 

0^  cos  (n,x)  + cos  (n,y)  * X (l*a) 

o COB  (n,y)  + T cos  (n,x)  ■ Y (l^b) 

y 


X and  Y ■ components  of  surface  forces  per  unit 
area  in  the  x and  y directions, 
respectively,  and 

cos  (n,x)  and  cos  (n,y)  ■ the  direction  cosines 


y 


b.  AN  ELEMENT  ON  THE  BOUNDARY 


Fig.  3.  Graphical  representation  of  equilibrium  and 
boundary  conditions 


Geometrically,  the  direction  cosines  can  be  represented  in  terms  of  dx  , 
dy  , and  ds  , as: 


cos  (n,y)  ■ - ^ 


l8.  Compatibility  condition.  A compatibility  condition  exists  if 
the  strain  components  at  any  point  A due  to  the  applied  loads  fit  their 
neighboring  points  perfectly  in  a manner  that  satisfies  strain  condition. 
The  compatibility  equation  in  terms  of  stress  components  for  plane 
strain  case  can  be  written  as 


3y7 


(o  0 ) * V (o  o ) 
X y X y 


V ■ Lapiaclan  operator 

19.  Stress  function.  As  stated  previously,  the  equilibrlujn  and 
compatibility  equation;,  represent  necessary  and  sufficient  conditions  to 


solve  for  thi-  stress  components  , and  . In  order  to 

solve  the  differential  equations  3 and  6,  a new  function  U(x,y)*  may  be 
introduced  such  as; 


In  order  to 


„ _ 3^0(x.y) 

°x  " 2 

3y 


xy  3x3y 


• This  function  U(x,y)  , which  is  called  stre"o  ’unction  or  Airy's 
stress  function,  was  first  introduced  in  1862  ty  G.  B.  Airy. 9 


V 


Adding  equations  Ta  and  Tb  yields: 


(8) 

^ 3y  3x 

20.  If  the  stress  components  , o^  , anJ  are  single- 

valued with  a continuous  first  and  second  derivative  throughout  the 
region  occupied  by  the  soil,  then  the  function  U(x,y)  must  have  a 
single  value  and  continuous  third-  and  fourth-order  derivatives  through- 
out the  soil  beneath  the  wheel. 

21.  It  is  apparent  that  the  equations  of  equilibrium  are  identi- 
cally satisfied  by  simply  substituting  equations  7a,  7b,  and  7c  into 
equations  3a  and  3b.  However,  a true  and  unique  solution  to  the  problem 
ceui  be  obtained  only  when  compatibility  is  also  satisfied.  This  condi- 
tion is  accomplished  by  substituting  equations  7a,  7b,  and  7c  into  equa- 
tion 6 as  follows: 


a‘*U(x,y)  , „ a‘*U(x.y)  , 3**U(x,y)  , ^ 

3x^  3x^3y^  3y^ 

Thus,  the  problem  of  solving  differential  equations  1,  2,  and  7 is  re- 
duced to  the  solution  of  one  blharmonlc  equation  that  satisfies  both  the 
equilibrium  and  compatibility  conditions.  The  solution  of  equation  11 
(given  later)  that  satisfies  the  boundary  condition  in  a simply  con- 
nected region*  is  a correct  solution. 


Complex  Representation  of  the  Weightless 
Plane  Strain  Problem 


22.  It  has  been  shown  that  the  evaluation  of  the  stresses  o , 

X 

a , and  t within  a simple  connected  region  lies  in  the  solution  of 

y 

a blharmonlc  differential  equation  such  as  equation  11  (given  below). 

The  solution  of  a biharmonic  equation  is  very  tedious  and  involved  with- 
out introducing  some  means  to  simplify  the  problem.  One  method  of 


* A region  R is  called  simply  connected  if  any  simple  closed  curve 
that  lies  in  R can  be  shrunk  to  a point  without  leaving  the  region. 
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Blmplifylng  the  problem  Is  to  transform  the  biharmonlc  equation  to  a 
harmonic  one*  using  some  analytic  function.**  Such  a procedure  Is  simi- 
lar to  that  developed  by  Muskhellshvlll''"*'  and  the  notations  used  are 

12 

those  adopted  by  Baladl . 

Representation  of  stress  function 
In  terms  of  harmonic  complex  function 


23.  Assuming  the  existence  of  a given  function  P( 
2 

lent  to  V U(x,y)  , equations  8 and  9 can  be  written  as: 

x,y)  equiva- 

0 ♦ 0 ■ V^U(x,y)  ■ P(x,y) 

X y 

(10) 

V^U(x,y)  ■ T'P{x,y)  ■ 0 

(11) 

Thus,  the  function  P(x,y)  Is  a harmonic  function  and  will  have  a con- 
jugate harmonic  function  Q(x,y)  . Consequently,  P(x,y)  + lQ(x,y)  Is 
the  analytic  function  of  z ■ x ♦ ly  and  can  be  written  as: 

W(z)  * P(x,y)  + lQ(x,y) 

(12) 

The  integration  of  this  function  with  respect  to  z yields  another  ana- 
lytic function.  Let  that  function  be  defined  as  U0[z)  . Consequently 

0(z)  ■ p(x,y)  ♦ iq(x,y)  ■ ^fw(z)  dz 

(13) 

where 

p(x,y)  and  q(x,y)  ■ real  and  Imaginary  parts,  respectively,  of 

0(z) 

Thus,  the  function  w(z)  can  be  obtained  by  differentiating  equation  13 
0'(z)  ■ i ■ F lQ(x,yj]  (!■*) 


* Harmonic  function  Is  the  one  that  satisfies  Laplace's  equation. 

**  Analytic  function  In  a region  Is  a function  whose  derivative  exists 
at  all  points  within  the  region. 
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2k,  Equating  the  real  parts  of  the  second  and  last  members  of 
equation  ll*  in  conjunction  with  the  Cauchy-Riemann*  condition  yields 


,15, 

or 

V^U(x.y)  - 2 - 2 - 0 (16) 

which  is  equivalent  to 

V^[u(x,y)  - xp(x,y)  - yq(x,y)]  *0  (17) 

25.  Equation  17  represents  a harmonic  function  equal  to  Pj^(x,y)  . 
Consequently,  the  stress  function  U(x,y)  is  equal  to 

U(x,y)  * xp(x,y)  + yq(x,y)  + Pj^(x,y)  (l8) 

Equation  l8  clearly  demonstrates  that  any  stress  function  can  be  formu- 
lated by  use  of  the  appropriate  conjugate  harmonic  functions  p(x,y)  , 
q(x,y)  , and  a harmonic  function  Pj(x,y)  • However,  the  functions 
p(x,y)  and  q(x,y)  were  used  for  a mathematical  expediency  but  not 
necessarily  for  the  .-.ohution.  This  necessitates  the  introduction  of 
another  harmonic  function  qj^(x,y)  conjugate  to  p^^(x,y)  , and  both 
functions  form  the  real  and  imaginary  parts,  respectively,  of  a function 
of  the  complex  variable  x(z)  where  z * x + iy  , such  that 

x(z)  = Pj^(x,y)  + lqj^(x,y)  (19) 

Thus,  the  stress  function  U(x,y)  can  be  expressed  as  follows: 

U(x,y)  » xp(x,y)  + yq(x,y)  + Rex(z)  (20) 

* The  Cauchy-Riemann  equations  that  satisfy  the  analytic  function  0(z) 
are  3p(x,y)/3x  = 3q(x,y)/3y  and  3p(x,y)/3y  ■ -[ 3q(x,y)/3x]  . 
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where 


Re  = real  part  of  a complex  function 
The  above  equation  can  be  rewritten  as: 

U(x,y)  * Re[(x  - ir)p(x,y)  + i(x  - iy)q(x,y)]  + Rex(z)  (2l) 

However,  x - iy  is  conjugate  of  z , which  may  be  denoted  by  ‘z  . 

Thus 

U(x,y)  = Re|7[p(x,y)  + iq(x,y)]  + x(z)f  (22) 

Recall  that  p(x,y)  + lq(x,y)  = 0(z)  . Consequently 

U(x,y)  = Re^z0(z)  + x(z)]  (23) 

26.  Through  use  of  the  property  of  complex  variable, equa- 
tion 23  can  be  transferred  to  the  following  form: 

2U(x,y)  = J0(z)  + z0(z)  + x(z)  + x(z)  (2l») 

Taking  the  partial  derivative  of  equation  2i*  with  respect  to  y , multi- 
plying it  by  i , and  then  adding  to  the  product  the  partial  derivative 
with  respect  to  x yields: 

Miizl  + i e 0(z)  + (25) 

However,  taking  the  second  portion  of  equation  2U  with  respect  to  x 
and  y and  then  adding  the  results  yields: 

V^U(x,y)  = ^0'(z)  + 0'  (z)]  ■ l*Re|0'(z)]  (26) 

Comparing  equation  26  with  equation  10,  then 

V^U(x,y)  = + Oy  = l4Re[0'(z)]  ^27) 
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Boundary  stress  in 
terms  of  complex  functions 

27 • So  far,  the  sum  of  stresses  a * a has  been  expressed  In 

X y 

terms  of  complex  potential;  thus,  the  next  logical  step  is  to  represent 
the  state  of  stress  at  a point  in  terms  of  complex  functions.  It  has 
been  shown  in  equation  7 that 


. a^U(x.y) 

. a^u(x.y) 
3x^ 

_ _ a^u(x.y) 


(7a  bis) 


(7b  bis) 


xy  3x3y 

Also,  geometrical  relationships  presented  in  fig.  3b  Indicate  that 


cos  (n,x)  * cos  (s,y)  * 


(7c  bis) 


cos  (n,y)  ■ -cos  (a,x)  * - 


(28a) 

(28b) 


28.  In  addition,  the  boundary  conditions  presented  previously  in 
equation  are: 


X » 0^  cos  (n,x)  + cos  (n,y) 


Y ■ o cos  (n,y)  + t cos  (n,x) 

y jQr 


Substituting  equations  7 and  28  into  equation  yields: 

,2, 

3y' 


X « ^ + a^UU.y)  ±_  r3U(x,y)] 

jy2  ds  3x3y  ds  ds  ^ 3y  J 

r3U(x,y)1 

ds  L 3x  J 


y . 3^0(x.y)  dx  3^U(x,y)  dy 

3x2 


(Ua  bis) 
(Ub  bis) 


3x3y  ds 


(29a) 


(29b) 
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Equations  29a  and  29b  can  be  combined  Into  one  equation  b>  multiplying 
the  second  equation  by  1 and  adding  the  first  one.  Then 


X + lY 


d ratJ(x.y)  _ au(x.y)] 
ds  L 3x  J 


or,  for  convenience,  this  equation  can  be  written  as 


X + lY  ■ -1  tt  r ♦ 1 


] 


(30) 


ds  I 3x 

Comparing  the  above  equation  with  equation  25,  it  can  be  written  that 


(31) 


X + lY 


-i  ^0(z)  + z0*(z)  + PTzT] 


(32) 


29.  The  value  of  ds  in  equation  32  corresponds  to  both  dx  and 
dy  as  indicated  previously  in  fig.  3b.  If  ds  traverses  along  the 
y axis,  then 


ds  * dy  = -i  dz 

Also 

X - 0 

X 

and 

y « T 

xy 

Substituting  these  values  into  equation  32  and  taking  the  derivative 
with  respect  to  z yields  the  following  equation; 

°x  * ^^xy  * (2)  ~ z0"(z)  - x"(z)  (33) 

If,  on  the  other  hand,  ds  traverses  only  along  the  x axis,  then: 


ds  * dx  ■ dz 


Also 
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and 


Substituting  into  equation  32  and  taking  the  derivative  with  respect  to 
z yields: 


a - iT  - 0'(z)  + 0'(z)  + zFTzT  + x"(z)  (3*») 

30,  Simpler  forms  can  be  obtained  by  adding  and  subtracting  equa- 
tions 33  and  3^: 


+ Oy  ■ 2[0'(z)  + 0'<z)]  * J4Re[0'(z)]  (35) 

and 

Oy  - 0^  - 2ix^y  - 2[z0"(z)  + x"(z)]  (36) 

Changing  i to  -i  on  both  sides  of  equation  36  yields  the  alternative 
form 


Oy  - 0^  ♦ 2iT^  ■ 2[z0''(z)  ♦ x"(z)] 


(37) 


For  convenience,  let  0'(z)  * ♦(z)  and  x"(z)  ■ <>(z)  ; thus,  equa- 
tions 35  and  37  can  be  presented  in  the  following  form: 

(38a) 
(38b) 


0^  + Oy  ■ 2|«(z)  ♦ ♦{z) J ■ >*Re[*(z)] 
Oy  - + 2ix^  ■ 2p*'(z)  + t|/(z)] 


Equation  38  can  be  used  to  determine  the  stress  components  in  terms  of 
complex  potentials  4(z)  and  V>(z)  . Thus,  if  definite  functions  for 
4(z)  and  ^(z)  are  chosen,  the  state  of  stress  from  equation  38  can  be 
obtained. 


2l» 


y 
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Introduction  to  Conformal  Mapping 


31.  In  previous  parts  of  this  report.  It  has  been  shown  that  the 
state  of  stress  at  a point  In  a region  can  be  defined  by  two  analytic 
functions  with  continuous  first  and  second  derivatives.  Thus,  the  solu- 
tion of  any  particular  problem  lies  In  obtaining  analytic  functions  that 
satisfy  the  stresses  at  any  point  within  the  region,  Including  the  bound- 
aries. However,  most  practical  problems.  Including  the  one  under  con- 
sideration, do  not  have  a simple  geometric  boundary.  Therefore,  the 
solution  to  the  problem  Is  to  transform  the  region  under  consideration 
with  Its  complicated  boundary  to  another  one  with  a smooth  boundary, 
such  that  each  point  In  the  original  region  corresponds  to  only  one  point 
In  the  new  region.  Once  the  solution  has  been  established  In  the  new  re- 
gion, It  can  be  carried  back  by  the  Inverse  transformation  to  the  region 
of  the  original  problem.  A transformation  that  possesses  the  property 

of  preserving  angles  of  Intersection  and  the  approximate  Image  of  small 
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shapes  Is  said  to  be  conformal  and  the  technique  used  Is  called  con- 
formal mapping. 

Schwarz-Chrlstoffel  transformation 

32.  According  to  Rlemann's  mapping  theory, an  analytic  function 
exists  that  will  map  any  pair  of  simply  connected  regions*  conformally 
onto  each  other.  To  clarify  this  statement,  let  In  fig.  l^a  be  a 
simple  closed  curve  In  the  z plane  (z  * x + ly)  forming  the  boundary  of 
a region  . Let  In  fig.  4b  be  another  simple  closed  curve  In 
the  t plane  (t  ■ r + Is)  forming  the  boundary  of  a region  Rg  . Then 
there  exists  a function  t ■ f(z)  analytic  In  R^^  that  maps  each  point 
In  the  Interior  of  onto  the  Interior  of  Cg  . However,  the  deter- 
mination of  a transformation  function  that  maps  conformally  a compli- 
cated region  onto  one  with  simple  geometry  Is  not  readily  available  and 
may  require  a series  of  auxiliary  functions  between  the  original  region 
and  the  region  In  which  the  stresses  at  any  point  can  be  defined  In 

* If  a simple  closed  curve  lying  In  a region  can  be  shrunk  to  a point 
without  leaving  the  region,  then  that  region  is  called  simply 
connected. 
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teras  of  analytic  functions.  The  transformation  functions  used  in 
this  study  were  obtained  from  Kober's  Dictionary  of  Conformal 

lli 

Representations. 

33.  One  method  of  special  interest  to  this  study  is  the  trans- 
formation of  a region  inside  a polygon,  whose  sides  never  cross  one 
another,  in  any  plane  onto  another  polygon  or  on  the  upper  half  of 
another  plane.  To  Illustrate  this,  let  w^  , shown  in 

fig.  5a  be  the  vertices  of  a polygon  in  the  w plane  (w  ■ u Iv)  that 

correspond  to  Interior  radian  angles,  , . . . . Also  let 

r,  , r„  , . . . r (r,  < r«  ...  < r ) be  the  points  on  the  real  axis  of 
1 c n 1 £ n 

the  t plane  (t  ■ r + is)  that  correspond  to  the  vertices  w^  , w^  , . . . 
w^  of  the  polygon,  respectively  (see  fig.  l<b). 

3l*.  A transformation  that  maps  the  interior  of  the  polygon  of  the 
w plane  onto  the  upper  half  of  the  t plane  and  the  boundary  of  the  poly- 
gon on  the  real  axis  of  the  t plane  is  given  by 


. where 

M and  N > complex  constants  that  control  the  size  and  position 
of  the  polygon 

Equation  39  is  called  the  Schwarz-Christoffel  transformation. 

35*  The  above  transformation  is  in  effect  mapping  a polygon  of 
the  w plane  onto  a similar  polygon  in  the  t plane  in  such  a way  that  the 
open  sides  of  the  polygon  extend  to  r * of  the  real  axis  of  the 
t plane.  It  may  also  be  said  that  the  region  bounded  by  the  opening  of 
the  w plane  and  a semicircle  with  infinite  radius  that  occupies  the  en- 
tire upper  half  of  the  t plane  comprise  the  interior  of  another  polygon 
in  the  t plane.  Therefore,  the  Schwarz-Christoffel  transformation  maps 
conformally  the  region  bounded  by  the  polygon  « '^2  » * • • ''n 
the  Interior  of  the  polygon  r^^  , r^  , ...  r^  and  the  semicircle  with 
infinite  radius  that  covers  the  entire  upper  half  of  the  t plane. 
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Transformation  of  Cartesian  coor- 
dlnates  to  curvilinear  coordinates 


36.  If  the  boundary  of  a region  Is  described  by  a set  of  orthog- 
onal curves,  a case  that  Is  encountered  In  many  practical  problems  and 
also  In  using  Schvarz-Chrlstoffel  transformation,  then  It  Is  more  con- 
venient to  use  curvilinear  coordinates  rather  theui  Cartesleui  coordi- 
nates. The  geometrical  transformation  between  the  Cartesian  and  curvi- 
linear coordinates  may  be  described  as  follows:  let  c In  fig.  6 be  a 
vector  originating  at  some  point  z ■ f(ty  , where  t ■ r ♦ Is  . Also, 
let  c^  and  c^  be  the  projection  of  c on  the  x euid  y axes  and  c^ 

and  c be  Its  projection  on  the  r and  s axes,  respectively.  According 
® 11 

to  Muskhellshvlll , the  transformation  equation  can  be  expressed  as: 

c + Ic  * (c  + Ic  )e”^®  (1*0) 

r s X y 


where 

6 > angle  between  c and  c measured  counterclockwise 
° X r 

•»16 

In  the  above  equation,  e"  Is  the  only  unknown  that  needs  to  be  de- 
termined. 

16 

37.  To  determine  e , let  the  point  z be  given  a displace- 
ment dz  in  the  direction  of  r . The  corresponding  point  in  the 
t plane  will  exhibit  a displacement  dr  > 0 in  the  direction  of  r . 
Since 


z ■ f(r  + Is)  ■ f(t) 


(1*1) 


then 


dz  ■ f'(t)  ds 


(1*2) 


By  definition 


dz  ■ |dz|  e^®  ■ f'(t)e^®  ds 


(1*3) 
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Fig.  6.  Transformation  of  Cartesian  coordinates  to 
curvilinear  coordinates 


1 


Equations  k2  and  )j3  can  be  combined  to  obtain 


f'(t)  ds  ■ lf'(t)|  |ds| 


(W) 


or 


e 


10 


(1*5) 


Therefore 


-10 


e 


(1*6) 


Dividing  equation  1*5  by  equation  U6  yields 


210  f'(t) 

e = 

f’(t) 


(1*7) 


Finally,  the  transformation  from  Cartesian  to  curvilinear  coordinates 
can  be  completed  by  substituting  equation  1*6  into  equation  1»3  to  obtain 


+ ic 


f(t) 

ffW 


(< 


c ) 

y 


(1*8) 


Representation  of  stress  com- 
ponents  In  curvilinear  coordinates 

38.  It  has  been  shown  in  equation  38  that  the  Cartesian  compo- 
nents of  stress  can  be  presented  in  terms  of  complex  potentials  ♦(z) 
and  4i(z)  . If  it  is  necessary  to  use  curvilinear  coordinates,  then  the 
complex  potential  can  be  expressed  as  a function  of  t , and  z can 
also  be  expressed  In  terras  of  t by  using  equation  1*1.  Ultimately,  it 

is  possible  to  express  the  Cartesian  stress  components  o , o , t 

X y xy 

at  a point  in  the  z plane  in  terras  of  Oj.  t ®g  » corresponding  to 

the  same  point  in  the  t plane.  In  order  to  do  that,  let  the  stresses 

0,0,1  be  the  normal  and  shear  stresses  acting  on  an  element, 
r s rs  Q 

as  shown  in  fig.  7,  in  the  t plane.  According  to  Timoshenko,  the  rela- 
tionship between  a , o , i and  o , o , t can  be  ex- 
pressed  as  follows: 
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I 


( 53a) 
(53b) 


+ 0^  > ^'[*(t)  + T(tf]  - j] 

For  convenience,  equation  53b  can  be  written  as 

, .21, 

°s  r rs  f'(t)  f (t)  J 

Thus,  if  4(t)  . .Kt)  . and  f(t)  are  known,  the  stresses  , 

and  T can  be  defined,  which  in  turn  allows  the  determination  of 


rs 

0,0,  and  t 
X ’ y 


xy 


3»» 
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PART  III:  ANALYTICAL  SOLUTION  OF  THE  PROBLEM 

Ul.  The  boundary  of  the  problem  under  consideration  (fig.  l)  con- 
sists of  a portion  of  a circular  arc  connected  with  two  straight  lines 
along  the  x axis  extending  to  infinity.  The  circular  arc  simulates  the 
area  of  contact  between  the  wheel  (l.e.,  rigid  wheel  or  pneumatic  tire 
with  inflation  pressure  equal  to  or  higher  than  the  average  contact 
pressure)  and  the  soft  soil  of  infinite  depth.  The  circular  arc  is 
acted  upon  by  a unllom  radial  stress  N and  a uniform  tangential 
stress  T , where  N and  T could  be  any  ct>mbinatiori  of  normal  and 
shear  stresses  developed  on  the  contact  area  between  the  moving  wheel 
and  the  soil.  In  short,  the  wheel-soil  interaction  problem  is  reduced 
to  the  evaluation  of  stressei.  within  a semi-infinite  region  bounded  by 
a circular  arc  and  two  straight  lines  extended  to  infinity. 

Geometric  Transformation 


U2,  Since  the  boundary  of  the  region  along  the  x axis  consists  of 
a circular  arc  and  two  straight  lines,  the  boundary  is  piecewise,  i.e,, 
the  derivative  at  the  point  of  contact  between  the  arc  and  the  straight 
lines  does  not  exist.  Therefore,  the  solution  lies  in  transforming  the 
region  with  its  piecewise  boundary  to  a semi-infinite  region  with  smooth 
boundaries  so  every  point  in  the  actual  region  will  correspond  to  one 
and  only  one  point  in  the  new  region.  This  will  be  accomplished  by  the 
following  transformation. 


Transformation  of 

half  plane  with  circular 

segment  removed  onto  infinite  strip 

1*3.  Since  direct  transformation  from  the  region  under  considera- 
tion onto  the  half-space  is  not  readily  available,  an  intermediate 
transformation  using  an  auxiliary  plane  will  be  used.  In  fig.  8,  let 
the  region  R^^  , in  the  z plane,  with  points  c,  d,  and  e having  their 
abscissas  at  x * +1  , x ■ 0 , and  x ■ -1  , respectively,  be  on  the 
region  R2  which  forms  an  infinite  strip  of  width  K in  the  upper 


I 
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Fig.  8.  Transformation  of  half  plane  with  circular  segment 
removed  on  infinite  strip 


half  of  the  w plane.  The  value  of  K is  directly  related  to  the  depth 
of  slnkage  5 by  the  following  equation 


6 ■ cot  K ■ R - Vr^  - 1 (5h) 

The  relations  between  6 , K , and  the  radius  of  the  wheel  R are  pre- 
sented in  fig.  9. 

ItU.  The  mapping  function  necessary  to  transform  the  region  , 
point  by  point,  onto  region  R2  can  be  expressed  as 

2w  _ ^ 

z ■ tanh  w « * - - — (55) 

e ♦ 1 

or 

w « tanh~^  z ■ r In  ^ (56) 

X * z 


Verification  of  the 
transformation  function 

U5.  Equation  55  can  be  easily  verified  by  resolving  z and  w 
to  their  real  and  imaginary  parts,  which  after  some  trigonometric  manip- 
ulation can  be  written  as 


X + iy 


tanh  u 


2 2 2 
cos  V tanh  u sin  v 


+ i 


cot  V sech  u 

2 2 
cot  V + tanh  u 


(57) 


Therefore,  by  substituting  the  values  of  u and  v at  any  known  point 
in  the  w plane,  the  corresponding  values  of  x and  y on  the  z plane 
can  be  obtained.  For  example,  at  point  c,  u » » and  y ■ 7t/2  j thus, 
the  corresponding  values  of  x and  y can  be  obtained  by  directly  sub- 
stituting the  values  of  u and  v in  equation  57.  Thus 


X + iy 


cos 


TfJ 


tanh  («») 

o p 

+ tan  (®)  sin 


(I) 


+ 1 


cot  sech^  (») 


cot 


(I) 


tanh“^  (■) 


(58) 
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Fig.  9.  Relationship  between  6 , R , and  K 
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y 


1 


1 


or 


X + ly  ■ 1 + i(0)  (59) 

I46.  By  equating  the  real  part  to  be  equal  to  the  real  part  and 
the  Imaginary  part  to  be  equal  to  the  imaginary  part  on  both  sides  of 
equation  51,  it  follows  that  x ■ 1 , and  y ■ 0 . The  calculated 
values  of  x and  y represent  exactly  the  real  and  imaginary  parts  of 
point  c on  the  z plane.  In  a similar  manner,  the  rest  of  the  points  can 
be  checked;  therefore  z ■ tanh  w is  truly  the  correct  transformation 
function. 

Transformation  of  the  infinite 
strip  onto  the  semi-infinite  plane 

1*7.  Since  the  transformation  of  the  original  region  was 

mapped  onto  the  infinite  strip  in  the  w plane,  the  problem  now  is 

to  map  the  region  bounded  b.y  the  infinite  strip  onto  the  upper  half  of 
another  plane,  such  as  t plane  (t  = r + is).  Theoretically,  the  bound- 
ary of  the  infinite  strip  shown  in  fig.  10  can  be  considered  as  a poly- 
gon with  its  vertices  at  infinity;  thus,  the  Schwarz-Cristoffel  equation 
can  be  applied.  The  Schwarz-Cristoffel  transformation  will  map  the  re- 
gion Rg  conformally  onto  the  region  , which  forms  the  entire  upper 
half  of  the  t plane.  Substituting  the  proper  values  of  a and  r of 

I 

points  a,  c,  d,  and  e obtained  from  fig*  10  into  equation  39  yields 


""L/"  1) " ” “ 2 + 0 


M and  N ^ complex  constants 

1»8.  To  determine  the  value  of  N , consider  the  condition  cor- 
responding to  point  a (i.e.,  t » <»  for  z * xi/2)  and  substitute  in 
equation  60  to  obtain  N » xi/2  . Also,  at  point  d,  t » 0 for 
z = i[{ir/2)  - K]  , which  after  substituting  in  equation  60  yields  a 
value  of  M = -2K/ii  • Once  the  values  of  M and  N are  known,  equa- 
tion 60  can  be  written  as: 


. y 


V 


} 


K , /t  - 1\  ^ Hi 

7 * - 


It  is  possible  to  eliminate  the  auxiliary  plane  w by  simply  combining 
equations  6l  and  thus 


Geometrically,  equation  62  means  that  the  boundary  of  the  z plane  with 
its  circular  indentation  will  be  mapped  on  the  r axis  of  the  t plane  and 
each  point  within  the  region  under  consideration,  , will  be  mapped 
on  a corresponding  point  in  the  upper  half  of  the  t plane. 

1»9.  Because  of  computer  software  limitations,  equation  62  may  be 
better  written  in  the  following  form: 

^2[-(K/ff)ln(t-l/t+l)-Hrl/2]  _ ^ 

* * 2[-(K/it)ln(t-l/t+l)+fll/2)  ^ 


Boundary  Conditions 


Boundary  conditions  in  the  z Plane 


30.  In  order  to  simplify  the  discussion,  it  is  more  convenient  to 
present  the  boundary  condition  equations  in  general  form  first  and  then 
modify  them  to  accommodate  the  specific  problem  of  this  study.  There- 
fore, let  ..he  normal  and  shear  stresses  acting  on  the  x axis  of  the 
z plane  be  defined  as  N and  T ; thus,  the  condition  of  equilibrium 

requires  that  N «uid  T be  equal  to  o and  x , respectively.  In 

y xy 


more  general  form,  N , T , a , and  x 

y xy 

equation  as  follows: 


can  be  combined  into  one 


N + IT 


0 + ix 
y xy 


0+0 
X J 


0 - o + 2ix 
y X XV 


51.  The  value  of  o + ix  can  be  obtained  easily  by  adding 
y xy 

equations  38a  and  38b.  Then 


y 


Cauchy  Integral  Formula 


Formula 

55.  In  the  following  section,  a brief  discussion  of  the  Cauchy 
Integral  formula  and  some  of  its  properties  relevant  to  this  study  is 
presented.  However,  more  detailed  study  on  this  subject  can  be  found  in 
any  standard  textbook  dealing  with  analytic  functions  of  complex 
varlables.^^*^^*^^  Letting  H{t)  be  some  analytic  function  anywhere 
within  the  boundary  C of  a simply  connected  region  R and  letting 
t^  be  any  point  within  the  region,  then  the  Cauchy  integral  can  be 
written  as 


1 

2iii 


/ 


(71) 


where  the  integral  Is  considered  positive  when  counterclockwise  around 


C . 


56.  The  Cauchy  integral  theorems  state  that  if  H(t)  is  a known 
function  on  the  boundary  of  a simple  closed  curve  C , then  the  value  of 
the  function  and  all  its  derivatives  can  be  found  at  any  point,  such  as 
t^  , inside  the  region  R ; then 


H(to) 


(72) 


As  a result  of  equation  72,  the  following  properties  of  Cauchy  integral 

, , 16 

are  in  orler: 

a.  If  the  function  H(t  ) is  harmonic  and  continuous  inside 

— o . 

the  region  and  ■ 0 , then  l/2itl  ^^H(t)/t  - tjdt 

is  equal  to  H(tg)  inside  the  region  and  to  zero  outside 
the  region.  This  property  is  usually  referred  to  as 
Poisson's  integral  formula  for  half-space. 10 

b.  If  the  function  H(t  ) is  harmonic  and  continuous  outside 

“ O - 

the  region  and  * 0 , then  l/2iri  - tjdt 

is  equal  to  zero  inside  the  region  and  -H(to)  outside 
the  region. 
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Application  for 
determining  stresses 

57.  Referring  to  equations  69  and  70,  it  can  be  stated  that  the 
functions  ♦(r)  , f(r)  [♦'  (r)/ f ' (r)J  , and  [f  * (r)/ f ' (r)]i|)(  r ) are  the 
boundary  values  of  the  harmonic  functions  in  the  lower  half  of  the 

t plcme  *(t)  , jV(t)/ f (t)^*'  (t)  , smd  (]f’  (t)/f ' (t)J^(t)  , respec- 
tively. Also,  the  functions  ♦ (r)  , f(r)^*'  (r)/f  (r^  , and 

♦nT[FTrT/f'(r)]  are  the  boundary  values  of  the  harmonic  functions 
in  the  upper  half  of  the  t plane  ♦(t ) , f(t)[^*(t)/f'  (t)]  and 

^[FTtT/f’(t)]  , respectively.  Since  all  these  functions  of  ,t 
ish  at  infinity,  then  the  Cauchy  integral  can  be  applied. 

58.  Starting  with  equation  70,  then 


van- 


«D  «» 

f HZL  ar  . r dr  . r 

2wi  J r - t 2*i  J r - t 2vi  J f'(rWr  - 


I 


TT 


dr 


^ f f <rj«(r)  ^ ^ f (H  - IT) 

2wi  J f'(r)(r  - t)  ^ 2vi  J r - t 


dr 


(73) 


Since  t is  a point  in  the  lower  half  of  the  t plane,  then  the  first, 
third,  and  fourth  integrals,  according  to  the  Cauchy  integral,  must  be 
equal  to  zero,  and  the  second  integral  is  equal  to  i(t)  ; thus,  equa- 
tion 73  can  be  reduced  to 


♦(t) 


1 f (N  - iT) 

2wi  I r - t 


dr 


(7Ua) 


However,  the  values  of  N and  T are  nonzero  only  at  the  interval 
-1  £ r ^ 1 ; therefore,  equation  7*»  can  be  further  reduced  to 


♦(t) 


1 f (M  - iT) 
2wi  / r - t 


dr 


(7Ub) 


which,  after  integration,  becomes 


♦(t) 


» - iT 
2«i 


-(m) 


(7Uc) 


from  which  ^'(t)  can  be  obtained  by  simply  differentiating  with  re« 
spect  to  t ; then 


♦ '(t) 


39<  In  a similar  manner,  the  value  of  i|i(r)  was  obtained  by 
operating  on  equation  69;  thus 

f N.t  IT  . _1_  /■  iilL  dr  + ^ r SEL  dr 

2irl  / r - t 2si  J r - t " 2itl  y r - t 

m 00 

+ _i_  f f{r)  *'{r)  f _rjj±_  »(r)  ^ gj 

‘’'I  FT7T 

Since  the  fvmctlons  f(r)  and  f'(r)  are  defined  inside  the  region, 
then  f(r)  and  f ' (r ) are  not  analytic  Inside  the  region  under  con- 
sideration. Thus,  the  Cauchy  Integral  cannot  be  applied  on  the  last  two 
terms  of  equation  76.  Therefore,  let  the  functions  G(t)  and  G'(t) 
be  two  analytic  fvmctlons  inside  the  region  so  that  at  the  boundary  they 
are  equal  to  f(r)  and  f ' (r)  , respectively.  Replacing  fTrT  with 
0(r)  and  f ' (r)  with  G'(r)  in  equation  76  yields 

♦1  ••  • 

f NjLil  dr  - ^ dr  f dr 

2si  jf  r - t 2si.  J r - t 2ifi  J r - t 

m 00 

_1_  r G(r|  ♦'jr)  ^ ^ _1^  r fl(r)  dr  (77) 

2wi  J G'(r)  (r  - t)  ^ 2wi  J G^)  (r  - t)  “ ^ ' 

Zm  Zw 


60.  According  to  the  Cauchy  integral  theorem,  the  third  term  in- 
volving ♦ (r)  is  equal  to  zero;  thus,  equation  77  after  integration  is 


reduced  to 


;N 


For  convenience,  the  value  of  ♦(t)  obtained  through  equation  7*»  can  be 
substituted  in  equation  78,  and  after  rearranging  terms,  the  value  of 
the  function  i|)(t)  can  be  expressed  as  follows: 


iK(t) 


0(t)*'(t) 


(79) 


where  ♦'(t)  and  f'(t)  are  known,  but  G(t)  and  0'(t)  are  yet  to 
be  determined. 


Determination  of  the  Function  G(t) 

61.  In  the  previous  section,  it  was  assumed  that  the  functions 

f (t)  and  0(t)  are  equal  at  the  boundary  of  the  region  under  consider- 
ation. Inside  this  region  only  G(t)  is  defined,  while  f (t)  can  be 
determined  outside  the  region.  The  graphical  relationship  between 
f(t)  , f(t)  , and  G(t)  is  represented  in  fig.  11. 

62.  In  order  to  determine  the  value  of  function  G(t)  , it  is 
necessary  to  transform  the  region  under  consideration  into  another  re- 
gion with  simple  geometry  such  as  a semi-infinite  plane.  This  kind  of 
transformation  cannot  be  performed  directly;  thus,  an  indirect  but 
simple  approach  may  be  employed  by  first  transferring  the  region  onto 
the  w plane  and  then  transferring  the  new  region  in  the  w plane  to  the 
upper  half  of  the  t plane. 

63.  Let  the  region  in  G(t)  be  defined  in  the  z plane  and  desig- 

nated by  Rj^  ; this  region  can  then  be  mapped  point  by  point  on  an  infi- 
nite strip  of  width  equal  to  + (it/2)  in  the  w plane  as  shown  in 

fig.  12.  The  mapping  of  region  Rj^  onto  R^  can  be  performed  easily 
using  the  same  transformation  function  described  in  equation  55  with  the 
exception  that  the  quantity  K is  replaced  by  . The  quantities  K 
and  can  be  expressed  by  the  following  equation 

- I - K (80) 

6*t.  The  second  step  is  to  transform  the  Infinite  region  , in 

t.6 
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the  w plane,  onto  a Beml-lnfinlte  region.  This  can  be  accomplished  by 
applying  the  Schwarz-Christoffel  transformation  to  map  conformally  the 
infinite  strip  , in  the  w plane,  onto  the  upper  half  of  the  t plane, 
as  shown  in  fig.  13.  The  procedure  used  is  exactly  the  same  as  that 
discussed  earlier  in  Part  III,  which  in  the  final  analysis  yields  the 
following  transformation  function; 


w ■ 


(81) 


65.  The;  w plane  can  be  eliminated  completely  by  simply  substi- 
tuting equation  8l  into  equation  55;  then 


2 * tanh 


(82) 


Although  equation  82  describes  the  geometry  of  the  region  under  con- 
sideration, it  does  not  map  the  entire  boundary  point  by  point.  However, 
if  a new  complex  function  1s  introduced  so  that  its  imaginary  part  is 
equal  to  the  imaginary  part  of  w while  its  real  part  is  equal  to  the 
real  part  of  w multiplied  by  K/(n  - K)  , then  the  region  Rj^  in  the 
z plane  will  be  mapped  point  by  point  on  the  t plane.  The  new  function 
Wj^  can  be  written  as  follows: 


+ ilm(w) 


(83) 


where  Im  is  the  imaginary  part  of  the  function.  If  the  function  w^^ 
is  substituted  into  equation  55,  the  value  of  G(t)  can  be  obtained  as 
follows : 


z * tanh  Wj^  = G(t)  (81») 

66.  From  equation  8U,  the  function  G’(t)  can  be  obtained.  Know- 
ing the  values  of  functions  G(t)  and  G'(t)  , the  function  i|i(t) 

(equation  79)  can  be  evaluated.  Finally,  the  stresses  0,0,  and 

X y 

T can  be  presented  in  the  following  equations: 
xy 

1*9 


s 


°x  " 


2Re^*(t)]  + Re^f{t) 

2Re[«(t)]  - RejTrrr 
T^y  - imjTTtT  ♦ <-(t)j 


Equations  85,  86,  and  8?  provide  the  solution  to  the  problem. 


Reduction  of  Data 


67.  The  computation  of 
the  WES  GE-UOO  computer. 


0 , and  T was  programmed  on 

y ^ 


Illustrative  Examples 


68.  Two  examples  are  presented  to  illustrate  the  use  of  the  theo- 
retical solution.  In  the  first  example,  a hypothetical  wheel  of  radius 
R is  partially  burled  in  soil,  as  shown  in  fig.  ll*,  so  that  the  effec- 
tive contact  length  and  the  sinkage  were  assumed  to  be  2 and  0.25  units, 
respectively.  The  wheel  was  also  assumed  to  exert  only  uniform  radial 
stress  along  the  contact  surface  between  the  soil  and  the  wheel.  A com- 
puter program  was  prepared  to  ceilculate  the  stresses  a , a , r , 

y x xy 

and  at  any  point  within  the  vicinity  of  the  wheel.  For  illustra- 

tion only,  the  stresses  along  a plane  1.2  units  below  the  soil  surface 
are  presented  in  fig.  ih  in  the  form  of  stress  distribution  diagrams  in 
the  interest  of  generality.  In  the  second  example,  the  same  hypothetical 
wheel  was  used  except  that  the  wheel  exerted  only  uniform  shear  stress 
along  the  contact  surface  between  the  wheel  and  the  soil.  The  stress 

distribution  diagrams  for  o , o , t , and  x are  presented  in 

y X mfix 

fig.  15. 


PART  IV:  LABORATORY  EXPERIMENTS 


69*  Ir.  previous  psrts  of  this  report,  a set  of  analytical  equa- 
tions we  developed  to  calculate  0,0,  and  t at  any  point 

X y xy 

within  the  soil  medium  in  the  vicinity  of  the  wheel.  These  stresses 
were  expressed  in  terms  of  the  radial  stress  N and  the  tangential 
stress  T that  might  he  generated  on  the  contact  surface  between  the 
soil  and  the  wheel  of  a moving  vehicle.  However,  since  the  radial  and 
tangential  stresses  are  greatly  Influenced  by  the  physical  properties  of 
both  wheel  and  soil,  the  study  of  the  stress  distribution  cannot  be  con- 
sidered complete  unless  it  is  associated  with  soil  strength  character- 
istics along  a potential  failure  surface  within  the  framework  of  the 
soil-wheel  system. 

70.  The  knowledge  of  the  physical  properties  of  soils  is  essen- 
tial not  only  to  estimate  the  relationship  between  the  normal  and  shear 
stresses  that  a soil  can  sustain,  but  also  to  check  the  capacity  of  the 
soil  to  support  the  vehicle  under  consideration.  It  might  be  worthwhile 
to  mention  that  vehicle  traction  Involves  the  topsoil,  which  does  not 
exceed  a few  inches  in  depth,  while  the  slnkage  problem  involves  the 
yielding  of  deeper  layers, 

71.  In  the  study  reported  herein,  a plastic  clay  was  tested  under 
controlled  conditions  to  simulate  field  conditions  as  closely  as  possi- 
ble. The  data  obtained  from  the  experimental  tests  were  then  applied  in 
the  theoretical  solution.  Because  the  objective  of  the  laboratory  tests 
was  to  simulate  field  soil-wheel  Interaction,  the  drainage  conditions 
and  the  speed  at  which  the  soil  beneath  the  wheel  is  deformed  were  also 
considered. 


Testing  Equipment  and  Material 


Equipment 

72.  The  entire  testing  program  was  performed  using  the  WES  ipl+7 
low-capacity  annular  shear  apparatus.  This  apparatus  is  described  in 
detail  in  reference  17.  The  test  specimen  is  U.50  in.  in  outside 
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diameter  and  2.30  in.  In  Inside  diameter,  giving  a cross-sectional  area 
of  11  sq  in.  (71  sq  cm).  Since  the  laboratory  testing  program  of  this 
study  was  intended  to  simulate  field  conditions  as  closely  as  possible, 
it  was  desirable  to  shear  the  soil  specimens  under  constant-volume  con- 
ditions at  a rate  of  strain  higher  than  that  commonly  used  in  soil 
testing. 

73.  Because  of  these  requirements,  the  torsion  shear  machine  was 
modified  and  adapted  to  conduct  torsion  shear  at  constant  volume  under 
vide  ranges  of  axial  stresses  and  rates  of  shear  deformation.  The  major 
modifications  that  were  introduced  were  the  constant-volume  loading  sys- 
tem and  the  automatic  recording  system.  Descriptions  of  the  constant- 
volume  loading  system  and  the  automatic  recording  system  are  presented 
in  Appendix  A. 

Material 

7^.  CH  soil  (Vicksburg  buckshot  clay)  was  used  in  this  testing 
program.  Index  properties  of  this  material  are  as  follows: 


Liquid  limit 

56 

Plastic  limit 

22 

Plasticity  index 

31* 

Activity 

0.88 

Specific  gravity 

2.68 

The  contents  of  one  carton  (80-85  lb)  were  thoroughly  mixed  on  a flat 
surface  and  then  split  into  5-lb  (mass)  portions,  which  were  stored  in 
plastic  bags  until  needed  for  tests. 

Sample  Preparation 

75.  The  soil  moisture-density  condition  used  in  the  testing  was 

X8 

determined  using  the  Modified  Berkeley  Pneumatic  Tamper.  Small 
batches  of  the  soil  were  mixed  at  different  water  contents  ranging  from 
21  to  28  percent.  The  material  was  mixed  by  adding  water  slowly  to  the 
dry  material  while  spreading  with  a spatula  to  break  down  any  large 
lumps  (this  material  was  not  forced  through  a screen).  The  soils  were 
stored  in  airtight  glass  Jars  and  allowed  to  cure  for  10  days.  After 

55 
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curing,  the  material  was  compacted  in  a standard  triaxial  mold  using  the 
tamper  set  at  11  lb  (force)  per  tamp,  with  30  tamps  per  layer  on  each  of 
eight  layers.  The  moistiure-denslty  curve  developed  for  these  compacting 
conditions  showed  an  optimum  water  of  25.6  percent  and  a maximum  dry 
density  of  95  pcf  (fig.  l6).  It  was  decided  to  prepare  specimens  for 
the  shear  tests  at  the  optimum  water  content  and  maximum  dry  density. 
Enough  material  for  35  shear  tests  was  mixed  at  the  proper  water  content 
in  a manner  very  similar  to  that  used  in  preparing  material  for  the  com- 
paction test. 

T6.  For  the  annular  shear  teats  on  soil  alone,  the  quantity  of 
material  required  to  produce  a compacted  specimen  3/l«  in.  thick  with  a 
dry  density  of  95  pcf  was  determined.  The  soil  was  placed  in  the  annu- 
lar shear  assembly  and  compacted  in  two  layers,  the  first  layer  being 
1/1»  in.  thick.  After  compaction  of  the  first  layer,  the  top  surface  of 
the  soil  specimen  was  scarified,  then  the  second  layer,  1/2  in.  thick, 
was  compacted. 

77.  The  pneumatic  tamper  was  used  with  a 3/8-ln.-dlam  Teflon 
tamping  foot  shaped  to  fit  a segment  of  the  annular  ring  (fig.  17).  The 
tamping  force  needed  to  produce  the  desired  specimen  thickness  was  cali- 
brated by  using  a number  of  tamping  forces  and  successively  correcting 
the  forces  until  the  desired  thickness  was  reached.  Three  coverages 
were  used  in  compacting  each  layer.  In  tests  to  determine  the  shearing 
resistance  of  soil  against  rubber  and  soil  against  steel,  one-half  of 
the  volume  of  a full  soil  specimen  was  used,  and  the  material  was  com- 
pacted in  a single  layer  with  three  coverages  of  the  tamper.  This  pro- 
cedure used  to  prepare  and  compact  the  soil  ensured  an  adequate  supply 
of  identical  specimens  for  the  testing  program,  and  it  is  believed  that, 
within  the  limits  of  experimental  variation,  the  specimens  tested  pos- 
sessed identical  initial  properties. 


Test  Program 


78.  Torsion  shear  tests  were  perfonned  on  four  types  of  speci- 
mens in  this  study,  as  follows; 
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Modified  Berkeley  pneiimatic  tamper 


a.  Horaogeneous . The  purpose  of  this  test  was  to  evaluate 
the  shear  strength  of  the  soil  Itself. 

b.  Rubber  on  soil.  In  this  type  of  test,  a smooth  rubber 
was  molded  in  the  form  of  a hollow  cylinder  to  occupy  the 
lower  half  of  the  torsion  shear  box,  while  the  upper  half 
was  filled  with  compacted  soil,  llie  purpose  of  this  test 
was  to  evaluate  the  shear  stress  developing  between  tire 
rubber  and  soil  under  different  radial  stresses.  The 
molding  of  the  rubber  ring  is  described  in  Appendix  B. 

£_.  Steel  on  soil.  This  test  was  similar  to  that  described 
in  subparagraph  b,  with  the  exception  that  the  molded 
rubber  was  replaced  by  polished  stainless  steel. 

d.  Monhomogeneous  soil.  In  this  test,  the  lower  part  of  the 
torsion  shear  box  was  filled  with  soil  compacted  at  a 
water  content  of  l6  percent,  the  surface  of  the  compacted 
soil  was  smoothed,  and  then  the  upper  part  of  the  box  was 
filled  with  the  soil  compacted  at  a water  content  of 
26  percent.  The  densities  of  the  upper  and  lower  layers 
after  compaction  were  about  95  and  89  pcf,  respectively. 
The  purpose  of  this  preliminary  test  was  to  initiate  a 
study  of  the  effect  of  differences  in  water  content  at 
the  shear  plane  on  soil  resistance. 

79.  The  initial  normal  stresses  used  in  each  type  of  test  were  5, 
15,  and  30  psl,  and  the  rates  of  shear  deformation  were  0.002,  0.2,  and 
2 in. /rain.  All  tests  were  conducted  under  constant-volume  conditions. 
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PART  V:  PRESENTATION  AND  DISCUSSION  OF  TEST  DATA 


80.  A summary  of  pertinent  Information  on  tests  performed  at  the 
time  of  this  reporting  is  presented  in  table  1,  and  the  shear  and  normal 
stresses  for  each  type  of  test  are  presented  graphically. 


Homogeneous  Soil  Specimens 


Test  results 

81.  The  CH  materials  compacted  in  the  annular  shear  box  at  an 
average  water  content  of  26  percent  and  dry  density  of  95  pcf,  grouped 
under  three  rates  of  shear  deformation,  are  presented  in  figs.  18-20. 

The  upper  curves  in  each  figure  represent  the  variation  of  the  shear 
stress  T with  respect  to  the  elapsed  time  during  the  test,  and  the 
lower  curves  represent  the  variation  of  the  normal  stress  o'  during 
shear  for  the  same  test.  The  elapsed  time  and  the  rate  of  shear  defor- 
mation in  each  test  were  used  to  obtain  data  on  the  amount  of  shear  dis- 
placement that  took  place  along  the  shear  plane. 

82.  Tests  for  which  results  are  shown  in  fig.  l8  were  conducted 
at  initial  normal  stresses  o^^  of  5,  15,  and  30  psi  and  sheared  at  a 
rate  of  shear  deformation  equal  to  0.002  in. /min.  The  same  initial 
normal  stresses  were  repeated  in  tests  whose  results  are  given  in 
figs.  19  and  20,  but  the  rates  of  shear  deformation  used  were  0.2  and 
2.0  in. /min,  respectively.  The  shear  stress,  fig.  l8,  increased  almost 
linearly  at  the  beginning  of  the  test  until  the  shear  deformation  was 
about  0.06  to  0.1  in.,  being  lower  for  tests  with  lower  initial  normal 
stress.  Then  the  shear  stress  decreased  gradually  until  the  average  de- 
formation was  about  0.2  in.  Beyond  that  point,  the  shear  stress  stayed 
constant  as  the  soil  behaved  as  a plastic  material.  The  pattern  of  nor- 
mal stress,  on  the  other  hand,  was  quite  different  from  that  of  shear 
stress.  Vfhlle  the  tests  conducted  under  of  3 psi  showed  an  in- 
crease in  the  normal  stress  during  shear,  indicating  the  tendency  of  the 
material  to  dilate,  the  other  two  tests  showed  a decrease  in  normal 
stress  with  increasing  deformation,  suggesting  a compressional  tendency 
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2.0  5.00  2.00  15  2.0  1*».55  2.97 
0.002  lU.OO  2.11  30  2.0  30.00  2.5U 
0.2  ll*.73  2.91 


Fig,  19.  Variation  of  shear  and  normal  stresses  for 
a homogeneous  soil  specimen  sheared  at  a rate  of  de- 
formation equal  to  0.2  in. /min 


Fig,  20.  Variation  of  shear  and  normal  stresses  for 
a homogeneous  soil  specimen  sheared  at  a rate  of  de- 
formation equal  to  2.0  in./min 


during  shear.  Consequently,  It  might  be  expected  that  negative  pore 
pressure  developed,  causing  an  Increase  in  strength  for  specimens  tested 
at  the  lower  Initial  normal  stresses,  while  at  higher  initial  normal 
stress.  It  might  be  expected  that  a positive  pore  pressure  developed, 
causing  a decrease  in  soil  strength  during  shear, 

83.  Because  of  the  nature  of  the  apparatus  used  In  this  study,  it 
was  not  possible  to  saturate  the  specimen.  It  Is  likely  that  the  pres- 
ence of  air  voids  In  the  compacted  specimen  made  any  changes  In  pore 
pressure  of  little  consequence.  Therefore,  assuming  no  development  of 
pore  pressures  from  the  beginning  of  the  test  up  to  the  failure  point,* 
total  stresses  would  be  effective  stresses.  Beyond  failure,  however, 

19 

the  shear  stresses  are  somewhat  affected  by  the  maximum  rate  of  flow, 
and  the  stress-displacement  curve  can  be  used  only  as  an  approximate 
curve. 

8t.  The  stress  versus  elapsed  time  curves  In  fig.  19  are  slightly 
different  from  those  presented  in  fig.  18  because  the  shear  deformation 
was  increased  from  0.002  to  0.2  In. /min.  As  shown  in  fig.  19,  shear 
stress  also  increased  linearly  until  the  shear  deformation  was  about 
0.1  in.j  then  it  started  to  curve  rapidly  until  it  reached  the  peak  and 
started  to  decrease  gradually  until  the  shear  deformation  was  about  3 in. 
Beyond  that  point,  the  shear  stress  showed  no  significant  change.  The 
curves  also  show  that  the  shear  deformations  at  failure  occurred  at 
about  0.6  to  0.8  in.,  being  lowest  for  tests  conducted  at  initial  normal 
stress  of  30  psi,  and  were  much  higher  than  those  shown  in  fig.  I8. 

Thus,  increasing  the  rate  of  shear  deformation  not  only  increases  the 
peak  shear  stress  but  also  the  shear  deformation  at  failure.  Increasing 
the  rate  of  shear  deformation  increases  the  tendency  of  the  soil  speci- 
men to  dilate,  which  is  reflected  by  the  great  increase  in  the  normal 
stress  for  specimens  sheared  at  initial  normal  stresses  of  5 and  15  psi. 
However,  the  test  conducted  at  an  initial  normal  stress  of  30  psi  indi- 
cated a tendency  toward  compression  during  shear. 


* The  state  of  failure  is  defined  by  the  point  of  maximum  shear  stress 
in  this  study. 
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83.  The  general  variations  of  shear  and  normal  stresses  with  re- 
spect to  the  elapsed  time  for  specimens  shown  in  fig.  20  are  similar  to 
those  presented  in  fig.  19 • although  the  peak  stresses  are  slightly 
higher,  reflecting  the  effect  of  increasing  the  rate  of  shear  deforma- 
tion. In  fig.  20,  the  shear  and  normal  stresses  for  specimens  sheared 
at  an  initial  normal  stress  of  30  psi  showed  a continuous  decline  beyond 
the  failure  point  simply  because  the  soil  started  to  extrude  from  the 
shear  box  and  the  assumption  of  constant  volume  was  not  valid  in  that 
particular  test.  The  same  thing  happened  when  the  test  was  repeated, 
"oil  strength  parameters 

86.  It  has  been  recognized  that  the  shear  strength  is  equal  to 
the  cohesion  c for  purely  cohesive  soil  and  is  equal  to  a tan  0 for 
purely  frictional  soil,  where  a is  the  applied  normal  stress  on  the 
failure  plane  and  0 is  the  angle  of  internal  friction.  However,  most 
soils  possess  the  properties  of  both  cohesion  and  friction;  for  these 
materials,  the  shear  strength  can  be  expressed  as  follows; 

T * c + o tan  0 (88) 

The  tests  described  so  far  in  this  report  to  simulate  in  situ  conditions 
were  of  relatively  short  duration  and  can  be  considered  unconsolidated- 
undralned  (l.e.,  Q tests)  for  all  practical  purposes.  These  tests  were 
conducted  to  serve  one  purpose:  to  obtain  a relationship  between  the 

shear  and  normal  stresses  that  may  occur  in  the  vicinity  of  a wheel  mov- 
ing on  a soil.  TT»e  strength  may  also  be  determined  directly  by  triaxial 
compression  tests,  and  results  can  be  correlated  with  cone  index,  which 
is  commonly  used  in  trafficablllty  studies.  However,  the  following  dis- 
cussion is  limited  to  compacted  soil  tested  with  annular  shear  apparatus 
and  sheared  at  a constant  rate  of  shear  deformation  under  constant 


87.  The  relationship  between  the  shear  stress  t and  the  normal 
stress  a for  the  CH  material  shown  in  fig.  21  siiggests  that  the  peak 
shear  stress  increases  with  increasing  initial  normal  stress  as  well  as 
the  rate  of  shear  deformation.  The  best-fit  line  for  the  experimental 
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Relationships  between  the  shear  and  normal  stresses  at  failure  for  homogeneous  soil  specimens 
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results,  as  obtained  by  the  method  of  least  squares,  yields  a cohesion 
intercept  equal  to  3.25  psi  and  (t  equal  to  31  deg.  By  substituting 
the  values  of  c and  0 in  equation  88,  the  linear  relationship  be- 
tween T and  0 can  be  expressed  as  follows: 

T - 3.25  + 0.620  (8S 


Initial 

Water 

Content.  % 

Dry 

Density 

pcf 

Cohesion 
c*  . psi 

Friction 
Angle 
0'  . deg 

Rate  of 
Strain 
in. /min 

26.0 

95.0 

3.25 

31 

0.002  to  2.0 

2l*.0 

98.0 

2.80 

21 

— 

2t.5 

93.1 

i*.17 

23 

0.03  to  0.6 

2I4.5 

92.3 

1.10 

29 

0.0036 

where  a is  measured  in  pounds  per  square  inch.  The  values  of  c and 
0 obtained  in  this  study  were  compared  with  those  obtained  on  the  same 
material  by  other  investigators  as  shown  in  the  tabulation  below: 


Type  of  Teat 


Q annular  shear  26.0  95.0  3.25  31  0.002  to  2.0 

on 

S direct  shear  2l*.0  98.0  2.80  21 

on 

R triaxial  2U.5  93.1  h.l?  23  0.03  to  0.6 

R trlaxial^^  2I4.5  92.3  1.10  29  0.0036 

88.  The  data  presented  in  the  table  indicate  that  the  value  c' 
obtained  in  this  study  is  within  the  range  of  those  found  by  other  in- 
vestigators; however,  the  value  of  0 was  slightly  higher  than  the 
values  obtained  by  direct  shear  or  triaxial  compression  tests.  This,  of 
course,  was  expected  because  the  average  rate  of  shear  strain  under 
which  the  tests  in  this  study  were  conducted  was  much  higher  than  the 
rates  used  in  other  studies.  The  difference  may  be  attributed  to  the 
nonuniform  strain  distribution  associated  with  the  annular  shear  test 
and  also  to  the  slight  change  in  pore  pressure  during  the  tests.  Some 
caution  must  be  applied  in  using  the  data  obtained  in  this  study  to  pre- 
dict the  strength  of  the  CH  material  as  accepted  by  standard  procedures. 
However,  these  data  provide  a more  realistic  picture  of  the  relationship 
between  x and  o in  soil  under  the  action  of  a moving  wheel  than  the 
same  data  obtained  from  other  tests. 

Effect  of  rate  of  shear 
deformation  on  peak  strength 

89.  Plastic  material  similar  to  viscous  material  exhibits  a 
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resistance  to  shear  strain  that  varies  with  the  speed  at  which  the  shear 
Ing  strain  is  applied.  For  perfectly  viscous  fluid  (l.e.,  Newtonian), 
the  applied  stress  is  directly  proportional  to  the  rate  of  shear  strain 
t and  can  be  expressed  as  follows: 

0 ■ (90) 


where 

n ■ the  viscosity 
£ ■ de/dt 

Since  clay  soil  is  a three-phase  system,  its  deformation  characteristics 
are  greatly  affected  by  environmental  conditions.  At  high  water  content 
clay  behaves  as  a viscous  fluid  similar  to  that  of  Newtonian  material; 
as  water  content  is  gradually  reduced,  the  flow  properties  of  clay  are 
also  reduced  and  its  plastic  properties  begin  to  disappear.  As  the 
water  content  decreases  further,  the  clay  soil  becomes  harder,  and  its 
properties  take  on  properties  of  solids  similar  to  those  of  Hookean  ma- 
terial where  time  has  no  effect  on  strength.  Under  most  natural  condi- 
tions, clay  soils  are  quite  different  from  Hookean  material  since  their 
strengths  depend  on  the  rate  of  strain  and  they  also  differ  from  Newto- 
nian material  because  their  viscosity  is  not  constant.  According  to 
19 

Hvorslev,  the  coefficient  of  viscosity  also  depends  upon  the  maximum 
velocity  gradient  and  upon  the  time  elapsed  alter  the  maximum  velocity 
has  been  attained. 

90.  The  relationship  between  the  shear  stress  t and  the  initial 
normal  stress  , as  shown  in  fig.  22,  suggests  that  the  strength  of 
the  CH  material  increases  with  Increasing  for  all  three  rates  of 

strain  used.  The  effect  of  the  rate  of  shear  deformation  on  the 
strength  of  the  soil  is  demonstrated  in  fig.  23,  which  clearly  shows 
that  there  is  a substantial  increase  in  the  strength  of  the  compacted 
CH  material  as  a result  of  increasing  the  rate  of  strain  from  0.002  to 
2.0  in. /min.  The  increase  in  strength  with  respect  to  I is  the  most 
significant  for  specimens  tested  at  an  initial  normal  stress  of  5 psl 
and  the  least  significant  for  those  tested  at  30  psl.  Such  an  increase 
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In  strength  Is  mostly  due  to  the  viscous  properties  of  clay  as  explained 
previously. 

91.  The  Increase  In  the  strength  of  soil  as  a result  of  Increas- 
ing the  rate  of  shear  deformation  has  been  reported  previously  by 
Taylor  (on  Boston  blue  clay)  and  also  by  Donaghe  (on  CH  material). 

In  both  Investigations,  trlaxlal  compression  apparatus  was  used;  however, 
the  trend  of  the  results  was  similar  to  that  presented  In  fig.  22.  As  a 
result  of  his  experimental  study,  Taylor  proposed  the  following 
expression: 

T ■ (o  + p^)|tan  0'  + (91) 

where 

Pj^  ■ the  Intrinsic  pressure  equal  to  c/tan  0 
f(de/dt)  ■ a function  of  £ 

Smooth  Rubber  on  Soil 


Test  results 

92.  In  this  series  of  tests,  the  CH  material  was  compacted  on 
smooth  tire  rubber  that  filled  the  lower  part  of  the  annular  shear  box. 
The  compaction  effort  was  similar  to  that  used  in  the  previous  series  of 
tests.  The  initial  normal  stresses  and  rates  of  shear  deformation  used 
In  this  series  were  exactly  the  same  as  those  for  the  homogeneous  soil 
specimens.  Results  of  the  constant-volume  tests  of  soil  shear  on  rubber 
In  the  annuleu*  shear  apparatus  are  presented  in  figs.  2b-26. 

93.  Fig.  2k  shows  the  variation  of  shear  and  normal  stresses  dur- 
ing the  tests  of  three  specimens  tested  at  a rate  of  shear  deformation 
of  0.002  In. /min  and  under  Initial  normal  stresses  of  1$,  and  30  psl. 
The  shear  stress  for  each  test  Increased  almost  linearly  with  respect  to 
time  until  failure.  Beyond  the  peak  point,  the  shear  stress  dropped 
slightly,  then  stayed  about  constant  till  the  end  of  the  test.  The 
normal  stress,  on  the  other  hand,  decreased  for  all  tests  with  increas- 
ing shear  deformation  from  the  beginning  of  the  test  until  slightly 
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after  the  peak  point;  then  it  leveled  off  at  about  the  same  time  that 
the  shear  stress  stayed  constant.  It  is  apparent  that  all  three  speci- 
mens, even  those  tested  at  an  initial  normal  stress  of  5 psi,  showed  a 
tendency  to  compress  or  consolidate  during  shear. 

9^.  Visual  examination  of  specimens  after  tests  showed  that  all 
the  movement  occurred  at  the  plane  of  contact  between  the  rubber  and  the 
soil,  which  was  characterized  by  a shiny  polished  surface  where  the  soil 
specimen  moved  as  an  intact  block  over  the  rubber.  Figs.  25  and  26  are 
graphs  similar  to  those  presented  in  fig.  2k,  with  the  exception  that 
rates  of  shear  deformation  used  were  0.2  and  2.0  in. /min,  respectively. 

A comparison  of  figs.  2k,  25,  and  26  shows  that  the  failure  point  became 
more  pronounced  with  increasing  rate  of  shear  deformation,  and  also  that 
the  shear  strengr.t.  Increased  with  Increasing  speed  of  the  test.  However, 
the  normal  stress  during  shear  did  not  seem  significantly  affected  by 
the  speed  of  the  test. 

Strength  parameters 

95.  Since  the  shear  plane  occurs  between  the  soil  and  the  rubber, 
the  concept  of  failure  based  on  Mohr-Coulomb  criteria  as  described  by 
equation  88  is  not  valid.  However,  it  might  be  possible  to  express  the 
shear  stress  occurring  at  the  plane  of  contact  between  the  soil  and  rub- 
ber by  the  adhesion  a and  surface  friction  angle  6 in  a manner  simi- 
lar to  that  of  Mohr-Coulomb  equation  for  homogeneous  material.  Then 

T ■ a + 0 tan  6 (92) 

96.  The  relationship  between  the  shear  and  normal  stresses  at 
failure  for  rubber  on  soil  is  presented  in  fig.  27.  The  figure  shows 
that  the  shear  stress  Increases  with  increasing  initial  normal  stress  as 
well  as  with  increasing  rate  of  shear  deformation.  Both  adhesion  and 
friction  Increased  with  increasing  rate  of  shear  strain.  The  average 
value  of  adhesion  increased  from  about  1.6  to  3.1  psi,  and  the  angle  6 
increased  from  12  to  19  deg  when  the  rate  of  strain  was  increased  from 
0.002  to  2.0  in. /min.  Therefore,  it  is  not  possible  to  correlate  the 
shear  strength  with  the  normal  stress  at  failure  in  a manner  similar  to 
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that  described  by  equation  92.  However,  if  the  best-fit  line  is  con- 
sidered in  expressing  the  relationship  between  t and  o at  the  inter- 
face between  soil  and  rubber,  then  the  following  relationship  can  be 
introduced: 


T * 2.3  + O.2650 


(93) 


Effect  of  rate  of  shear 
deformation  on  strength  at 
interface  of  soil  and  rubber 

97.  The  relationship  between  the  maximum  shear  stress  developed 
at  the  interface  of  the  rubber  and  the  CH  material  and  the  rate  of  shear 
deformation  £ is  depicted  in  fig.  28.  There  is  a striking  resemblance 
between  this  figure  and  fig.  23,  both  showing  that  the  shear  stress  in- 
creased with  increasing  rate  of  shear  deform-tlon.  Both  figures  show 
that  the  curve  for  each  set  of  tests  Is  concave  upward,  giving  the  im- 
pression that  the  shear  stress  at  failure  Increases  indefinitely  with 
respect  to  the  increase  of  the  rate  of  strain.  However,  it  is  also  pos- 
sible that  the  curve  might  change  course  and  follow  another  direction 
with  increasing  rates  of  strain  beyond  those  used  in  this  test  program. 
This  point  needs  further  study. 

Effect  of  wetness  on 
the  shear  stress  at  the 
Interface  of  soil  and  rubber 

98.  Three  tests  were  conducted  in  a manner  similar  to  that  used 
in  previous  tests  with  the  exception  that  the  rubber  surface  was  covered 
with  a film  of  water,  using  wet  paper  toweling  prior  to  preparing  the 
specimen,  and  all  specimens  were  sheared  at  a rate  of  deformation  of 
2.0  in. /min. 

99.  The  values  of  the  shear  and  normal  stresses  at  failure  of  the 
wetted  rubber  surface  are  presented  in  fig.  29;  on  the  same  figure  re- 
sults of  comparable  tests  with  a dry  inibber  surface  are  also  shown. 

Fig.  29  clearly  shows  that  wetting  the  rubber  surface  reduced  the  shear 
resistance  to  about  50  percent  of  what  it  was  before  wetting,  even 
though  the  normal  stress  was  not  appreciably  affected.  This  may  explain 
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why  vehicles  with  pneumatic  tires  become  immobilized,  even  after  rain  of 
short  duration. 


Steel  on  Soil 


Test  results 

100.  In  this  teat  series,  the  CH  material  obtained  from  the  same 
batch  used  in  previous  tests  was  crn^pacted  on  a smooth  and  polished 
stainless  steel  that  filled  th.-  icwer  part  of  the  annular  shear  box. 

The  test  was  conducted  under  thi;  stuae  conditions  that  were  imposed  dur- 
ing the  rubber  on  soil  series. 

101.  Plots  of  the  shear  and  normal  stresses  versus  elapsed  time 
with  test  data  grouped  under  each  of  the  three  rates  of  shear  deforma- 
tion used  in  the  test  program  are  shown  in  figs.  30-32.  The  variations 
in  shear  and  normal  stresses  with  respect  to  time  of  these  three  groups 
greatly  resemble  those  presented  in  figs.  2U-26  for  the  rubber  on  soil 
series.  The  average  deformation  to  failure  was  about  the  seime  for  both 
test  series,  while  the  peak  shear  stress  was  slightly  lower  for  steel  on 
soil  than  for  rubber  on  soil. 

Strength  parameters 

102.  The  relationship  between  the  shear  and  normal  stresses  at 
failure  for  soil  specimens  sheared  on  the  polished  stainless  steel  is 
shown  in  fig.  33.  The  figure  shows  that  the  peak  shear  stress  Increased 
with  Increasing  initial  normal  stress  and  also  with  the  speed  at  which 
the  test  was  conducted.  However,  the  spread  of  the  results  is  much  less 
than  that  observed  for  the  rubber  on  soil  series  shown  in  fig.  27.  Also, 
the  strength  of  soil  sheared  on  polished  steel  was  about  half  the 
strength  of  soil  sheared  on  rubber  under  comparable  testing  conditions. 
Therefore,  it  may  be  concluded  that  the  traction  forces  associated  with 

a pneumatic  tire  are  larger  than  those  associated  with  a rigid  wheel. 

103.  The  best-fit  line  of  the  failure  points  shown  in  fig.  33  in- 
dicates an  adhesion  intercept  of  1.1  psi  and  an  angle  of  surface  fric- 
tion equal  to  9.5  deg.  Based  on  the  experimental  values  of  a and  6 
obtained,  the  linear  relationship  between  the  shear  and  normal  stress  at 
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CLAFSCO  TIME,  MIN 

b.  NORMAL  STRESS 


Fig.  30.  Variation  of  t and  o with  time  for  soil 
sheared  against  polished  steel  at  rate  of  shear  de- 
formation of  0.002  in. /min 


ELAPSED  TIME,  MIN 


b.  NORMAL  STRESS 

Fig.  31.  Variation  of  x and  o with  time  for  soil 
on  polished  steel  at  rate  of  shear  deformation  of 
0.2  in. /min 
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Relationship  between  normal  and  shear  stresses  at  failure  at  interface  of  soil  and  steel 


the  contact  surface  between  soil  and  polished  steel  can  be  expressed  as 
follows : 


T ■ 1.1  + 0.15c 


(91*) 


Effect  of  rate  of  shear 
deformation  on  strength  at 
interface  of  soil  and  steel 

lOl*.  The  relationships  between  the  maximum  shear  stress  that  de- 
veloped at  the  interface  of  the  soil  and  the  polished  steel  and  the  rate 
of  shear  deformation  are  shown  in  fig.  3l*.  Once  again  the  peak  shear 
stress  showed  an  increase  with  increasing  rate  of  shear  deformation. 
However,  the  rate  of  the  increase  was  much  lower  than  that  observed  at 
the  interface  of  the  soil  and  rubber  or  in  the  homogeneous  soil. 

Effect  of  wetness 
on  the  shear  stress  at 
interface  of  soil  and  steel 

105.  Three  tests  were  conducted  in  which  the  surface  of  the  pol- 
ished stainless  steel  in  contact  with  soil  was  covered  with  a film  of 
water  in  a manner  similar  to  that  used  in  the  rubber  on  soil  series. 

The  tests  were  also  conducted  at  initial  normal  stresses  of  5>  15,  and 
30  psi,  and  all  were  sheared  undrained  at  a rate  of  shear  deformation 
equal  to  2.0  in. /min. 

106.  The  relationship  between  the  shear  and  normal  stresses  for 
specimens  tested  with  a wet  steel  surface  and  those  with  a dry  steel 
surface  is  shown  in  fig.  35.  The  figure  clearly  shows  that  the  water 
film  decreased  the  peak  shear  stress  to  about  65  percent  of  that  with  a 
dry  steel  surface.  Thus,  the  degree  of  wetness  would  appreciably  reduce 
the  traction  characteristics  of  a rigid  wheel. 

Nonhomogeneous  Soil  Specimens 


Test  results 

107.  The  CH  material  in  this  test  series  was  compacted  in  two 
layers  of  equal  thickness  but  with  two  different  degrees  of  wetness. 
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and  steel 


The  water  contento  ol’  the  upper  and  lower  layers  were  26  and  l6  percent, 
respectively,  and  both  layers  were  compacted  with  the  same  compaction 
effort  used  in  the  homogeneous  soil  series.  Other  testing  conditions 
such  as  the  initial  normal  stresses  and  the  rates  of  shear  deformation 
used  in  this  test  were  kept  the  same  as  those  used  in  previous  test 
series. 

108.  It  was  expected  that  the  stress  versus  deformation  curves 
for  the  nonhomogeneous  specimens  would  be  comparable  to  those  of  the 
homogeneous  specimens  on  the  assumption  that  the  portion  of  the  soil 
specimen  with  the  lower  water  content  would  have  the  major  effect  on 

the  strength  and  deformation  characteristics  of  the  nonhomogeneous  speci- 
mens. However,  the  results  showed  that  although  the  shape  of  the  curves 
(not  shown)  resembled  very  much  those  presented  in  figs.  l8-20,  the  peak 
strengths  of  the  nonhomogeneous  soil  specimens  were  much  higher  than 
those  determined  by  comparable  tests  of  the  homogeneous  soil  specimens. 
Soil  strength  parameters 

109.  The  relationship  between  the  shear  and  normal  stresses  at 
failure,  shown  in  fig.  36,  suggested  a cohesion  Intercept  of  ?.5  psi  and 
an  angle  of  internal  friction  of  37  deg.  It  is  apparent  that  the  value 
of  0'  based  on  the  best-fit  line  is  about  6 deg  higher  than  that  for 
the  homogeneous  soil,  while  the  cohesion  intercept  is  about  0.7  psi 
lower.  This  difference  may  be  attributed  to  the  slight  migration  of 
water  from  the  upper  layer  to  the  drier  lower  layer,  making  the  overall 
water  content  at  the  failure  plane  of  the  nonhomogeneous  specimen  lower 
than  26  percent,  which  caused  the  nonhomogeneous  soil  specimen  to  be- 
have in  a less  plastic  manner  than  the  homogeneous  specimens. 

110.  Because  of  funding  limitations,  it  was  not  possible  to  test 
nonhomogeneous  specimens  of  which  both  halves  were  quite  wet  of  optimum, 
with  one  half  very  much  wetter. 

Relationship  Between  the  Analytical  and 
Experimental  Investigations 

111.  In  Part  III  of  this  report,  it  is  shown  that  the  state  of 
stress  at  any  point  within  a soil  under  a uniformly  loaded  circular 
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segment  can  be  expressed  as  follows; 


Oy  ■ 2Re[*(t)J  + Re^f(t)  + i|;(t)j  (85  bis) 

0^  = 2Re[i|i(t)]  - Re^f(t)  + i|i(t)j  (86  bis) 

(87  bis) 


where  4>(t)  , f(t)  , and  i|/(t)  are  analytic  functions  that  depend  on 
the  geometrical  configuration,  boundary  condition,  and  applied  stresses. 
The  function  ij((t)  depends  on  the  values  of  4i(t)  and  f(t)  , which 
can  be  expressed  as  follows : 


. / . X N - IT  , 
»(t)  = -7rr. — In 


2x1 


(Hi) 


f(t)  = tanh 


k , /t  - l\  . xi 

- ¥ ^^\rTT)  iy 


(7**c  bis) 

(95) 


It  is  quite  obvious  that  the  function  f(t)  depends  only  on  the  geomet- 
rical shape  of  the  medium  (i.e,,  soil-wheel  system),  while  *(t)  de- 
pends on  the  applied  radial  and  tangential  stresses  N*  and  T ,*  re- 
spectively. 

112.  In  Part  IV  of  this  report.  It  la  shown  that  the  value  of 
shear  stress  t , which  Is  equivalent  to  T , can  be  linearly  related  to 
the  applied  normal  stress;  such  a linear  relationship  Is  greatly  influ- 
enced by  the  type  of  wheel  and  soil  properties.  Therefore,  if  the  type 
of  wheel  used,  the  properties  of  the  soil,  and  the  applied  normal 
stresses  are  known,  it  is  possible  to  estimate  the  value  of  the  corre- 
sponding shear  stress.  If  the  normal  and  shear  stresses  at  the  soil- 
wheel  interface  and  also  the  geometry  and  boundary  conditions  of  the 
soil-wheel  systems  are  known,  it  is  possible  to  evaluate  the  functions 


• The  values  of  N and  T in  the  analytical  solution  are  equivalent 
to  0 and  t , respectively. 
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4(t)  , f(t)  , and  i|i(t)  . Once  these  functions -are  known,  equa- 
tions 85-87  can  be  used  to  evalxiate  the  state  of  stress  at  any  point 
beneath  the  wheel.  The  procedure  outlined  is  Illustrated  by  the  fol- 
lowing example.* 

113.  Consider  a 51-ln.-diam  rubber-rimmed  wheel  Imbedded  3 in.  in 
soil  with  properties  similar  to  the  CH  material  used  in  this  study  (see 
fig.  37a).  Assuming  that  the  wheel  is  in  a state  of  impending  motion 
while  exerting  a uniform  radial  strers  of  25  psi  on  the  underlying  soil, 
it  is  desired  to  evaluate  the  state  of  stress  at  point  A 1.2  ft  below 
"the  soil  surface  and  2.3  ft  to  the  left  of  the  center  line  of  the  wheel. 

llU.  The  first  step  is  to  estimate  the  maximum  shear  stress  that 
can  be  developed  along  the  Interface  of  the  soil  and  the  wheel,  using 
equation  93  (t  * 2.3  + O.2650).  Then: 

T » 2.3  + (0.265  X 25)  or  8.925  psl 

115.  The  second  step  is  to  idealize  the  wheel  so  that  the  analyt- 
ical solution  can  be  applied  directly.  This  can  be  done  by  introducing 
such  a scale  factor  that  the  projection  of  the  contact  surface  L from 
the  center  line  of  the  wheel  is  equal  to  unity,  as  shown  in  fig.  37b. 

116.  The  third  step  is  to  evaluate  the  value  of  K using  either 
equation  5**  or  fig.  9-  Then 

^ = cot  K » 0.25 

K = 1.32587  radians 

Once  the  values  of  K , c , and  t are  known,  it  is  possible  to  evalu- 
ate the  stresses  directly  using  the  computer  program  in  Appendix  A to  get 
the  values  of  o , o , and  t for  the  point  under  consideration. 

For  this  problem,  the  computer  program  yielded  the  following  stresses  at 

point  A:  o = -2.6U6  pal;  o = -O.165  psi;  and  t = I.8II  psi. 

X y xy 

* The  example  is  presented  to  Illustrate  the  method  of  solution  pro- 

posed in  this  report  and  is  not  necessarily  the  solution  of  a practi- 
cal problem. 
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117 • Because  of  the  similarity  between  the  Ideal  wheel  presented 
In  fig.  37b  and  that  shown  In  figs.  ll<  and  13 • diagrams  presented  In 
figs.  1I4  and  13  can  be  used  to  evaluate  the  state  of  stress  at  point  A. 
This  can  be  done  simply  by  measuring  the  ordinates  of  respective  stress 
diagrams  and  multiplying  them  by  23  for  evaluating  stresses  resulting 
from  the  normal  stress  alone  and  by  8.923  for  evaluating  stresses  re- 
sulting from  shear  stress  alone.  The  resulting  stresses  can  be  added 
together  by  applying  the  principle  of  superposition  as  shown  below: 

a.  Stresses  at  point  A due  to  normal  stress  alone 

- -0.239  X 23  - -3.973  psi 

a ■ -0.0363  X 23  ■ -1.1»13  psi 

y 

T * O.12U  X 25  - 3.100  psi 
xy  ' e 

b.  Stresses  at  point  A due  to  shear  stress  alone 

o^  « 0.372  X 8.923  ■ 3.320  psi 

a « 0.1402  X 8.923  « 1.231  psi 

T * -O.lUU  X 8.923  * -1.283  psi 

xy 

By  adding  the  stresses  in  a and  b,  the  following  results  are  obtained: 

®x  " "5*975  + 3.320  • -2.633  psi 

o ■ -1.413  ♦ 1.231  “ -0.162  psi 

y 

T ■ 3.100  - 1.283  ■ 1.813  psi 
xy 

These  values  agree  closely  with  those  yielded  by  the  computer  program 
and  presented  in  paragraph  II6, 
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PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


ConcluBlons 

118.  The  following  eoncluBlona  are  believed  warranted  by  the 
findings  of  this  study; 

a.  The  state  of  stress  and  the  orientation  of  stresses  at 
any  point  within  a region  beneath  a wheel  can  be  ob- 
tained analytically  in  terms  of  the  radial  and  tangen- 
tial stresses  generated  by  the  wheel.  For  plane  strain 
cases,  the  normal  and  shear  stresses  beneath  the  wheel 
are  independent  of  the  properties  of  the  supporting 
medium. 

b.  The  soil-wheel  Interaction  problem  can  be  solved  by  a 
rational  approach  using  the  instantaneous  stresses  in 
providing  design  criteria  regarding  trafficability 
problems . 

c_.  The  experimental  study  demonstrates  that  much  pertinent 
Information  regarding  Immobilization  due  to  inadequate 
traction  can  be  obtained  through  relatively  simple  lab- 
oratory tests. 

d.  Based  on  limited  experimental  tests,  it  was  found  that 
the  traction  forces  on  any  potential  failure  plane  in- 
creased with  increasing  rate  of  shear  deformation  and 
increasing  normal  stresses  acting  on  that  plane.  How- 
ever, these  traction  forces  decreased  sharply  with  in- 
creasing wetness  of  the  potential  failure  surface. 

£.  The  maximum  shear  stress  at  the  Interface  of  soil  and 
smooth  rubber  is  higher  than  that  at  the  interface  of 
soil  and  polished  stainless  steel.  Stresses  at  both 
these  Interfaces  are  lower  than  the  strength  of  the  soil 
alone. 

f . The  relationship  between  the  shear  and  normal  stresses 
obtained  through  the  experimental  tests  can  be  combined 
with  the  analytical  solution  to  provide  rational  solu- 
tions in  the  field  of  terramechanics. 


Recommendations 


119.  It  is  known  that  the  study  of  relationship  between  soils  and 
wheels  of  moving  vehicles  could  take  various  approaches  depending  on  the 
type  of  vehicle  and  the  performance  desired.  With  regard  to  further 


development  of  the  findings  of  this  study,  the  following  recommendations 
are  made; 

a.  Extend  the  theoretical  solution  to  Include  cases  In 
which  the  radial  and  tangential  stresses  between  the 
soil  and  the  wheel  vary  along  the  surface  of  contact. 

The  variation  of  these  stresses  could  be  assumed  to  be 
sinusoidal,  parabolic,  or  other  geometric  distributions 
that  more  realistically  simulate  field  conditions. 

b.  Present  the  analytical  solution  of  wheels  with  different 
sizes  and  slnkage  In  the  form  of  Influence  charts  and 
tables,  thus  eliminating  the  need  for  tedious  calcula- 
tion or  Individual  computer  runs. 

c_.  Extend  the  present  solution  to  develop  expressions  for 
the  evaluation  of  strains  and  displacement  within  soil 
In  the  vicinity  of  the  wheel. 

d.  Follow  the  same  procedures  In  evaluating  stresses,  but 
assume  the  soil  to  be  a linear  viscoelastic  material. 

£.  Extend  the  scope  of  the  experimental  study  on  the  plas- 
tic clay  used  In  this  study  and  also  study  other  soils, 
particularly  with  respect  to  loss  of  traction  due  to 
sllpperlness. 

f.  Verify,  if  possible,  the  findings  of  this  study  using 
true  model  testing. 

Adapt  the  analytical  solution  for  the  study  of  stresses 
around  circular  cavities  at  or  near  the  soil  surface. 


120.  Other  recommendations  that  are  related  to  the  vehicle  trac- 
tion problem  but  are  not  directly  related  to  the  findings  presented  in 
this  study  are: 

a.  Study  the  plastic  equilibrium  of  soil  in  the  vicinity  of 
a wheel  under  the  Influence  of  forces  generated  by  a 
moving  vehicle.  This  problem  could  be  treated  by  either 
considering  the  distribution  of  normal  stresses  along 
the  rupture  surface  (i.e.,  stability  problem)  or  by  as- 
suming that  every  point  in  the  rupture  zone  Is  in  a 
state  of  limiting  equilibrium  (i.e.,  limiting  equilib- 
rium problem).  In  either  case,  use  must  be  made  of  an 
appropriate  failure  criterion  such  as  Mohr,  Tresca, 

Von  Miseo,  or  others. 

b.  Apply  the  procedure  outlined  in  subparagraph  a to  solve 
the  problem  of  a grouser  analytically  or  experimentally 
or  both,  taking  Into  consideration  the  effect  of  size 
and  shape  of  the  grouser  and  of  the  distance  between 
two  adjacent  grousers  on  the  performance  of  a vehicle. 
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Develop  as  an  ultimate  goal  a prediction  equation  that 
can  relate  the  energy  or  force  input  by  the  vehicle  to 
the  erergles  or  forces  lost  by  traction,  friction,  and 
slippac-e;  soil  compaction;  and  soil  flov  and  other  actual 
measureuents  or  reasonable  estimates  of  these  quantities. 
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APPENDIX  A:  ANNULAR  TORSION  SHEAR  APPARATUS 


1.  The  test  apparatus  used  in  this  program  was  an  annular  shear 
device  designed  and  built  by  the  U.  S.  Army  Engineer  Waterways  Experi- 
ment Station  (WES)  in  19^7.  The  device  is  described  in  detail  in  ref- 
erence 23.*  The  annular  shear  test  was  employed  in  this  Investigation 
because  any  desired  amoimt  of  displacement  may  be  obtained  with  an  an- 
nular test  specimen,  and  shear  strengths  at  large  displacements  were  of 
Interest.  The  19**7  version  of  the  annular  shear  apparatus  had  to  be 
modified  so  that  data  could  be  recorded  electronically;  this  modifica- 
tion was  necessary  since  the  rapid  rates  of  displacement  used  would  not 
permit  manual  observation. 


Shear  Assembly 

2.  The  unmodified  device  employed  a load  frame  to  measure  the 
torque  from  which  the  shear  stress  was  computed.  A section  through  the 
unmodified  shear  apparatus  is  shown  in  fig.  Al.  The  soil  specimen  is 
confined  by  inner  and  outer  rings  that  are  divided  into  an  upper  and 
lower  pair.  The  upper  inner  ring  is  guided  by  a special  ball  bearing, 
permitting  both  rotation  and  vertical  movement.  The  lower  ring  fits 
snugly  over  a rib  in  the  rotating  base,  which  is  supported  and  guided 

by  bearings  on  the  stationary  base.  The  rotating  base  can  be  rotated  by 
hand  or  motor  through  the  self-locking  worm  gear  that  holds  the  base 
stationary  when  the  drive  shaft  is  not  rotated. 

3.  For  this  testing  program,  the  shear  apparatus  was  modified  by 
replacing  the  upper  and  lower  dentated  porous  plates  with  nonporous 
plates  having  blades  0.03^  in.  thick  and  1/6  in.  high  and  spaced  at 
30-deg  Intervals  around  the  annulus.  The  bladed  plates  are  shown  in 
fig.  A2;  the  top  bladed  plate  is  shown  attached  to  the  torque  wheel. 

The  device  was  modified  in  this  manner  because  for  soil  on  soil  tests. 


• Reference  numbers  refer  to  like-numbered  items  in  the  list  of  refer- 
ences at  the  end  of  the  main  text. 
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annular  shear  apparatus 


the  dentated  porous  plates  could  not  provide  '.ufficient  frictional  re- 
sistance against  the  top  and  bottom  surfaces  of  the  soil  specimen  to 
prevent  slippage.  It  was  not  necessary  that  the  bladed  plates  be  porous 
since  undrained  tests  were  performed.  The  zones  of  disturbance  at  the 
top  and  bottom  of  the  soil  specimen  were  not  thought  to  influence  the 
test  results  since  the  blades  were  thin  and  short  (1/8  in.  high)  and  the 
soil  specimen  was  3/U  in.  thick. 

Torque  Loading  System 

4.  As  the  rotary  base  rotates,  a torque  is  produced  by  pulling  on 
two  flexible  wire  cables  that  are  attached  to  the  torque  wheel.  The 
cables  pass  over  horizontal  and  vertical  guide  sheaves  and  connect  to 
the  balance  beam,  which  reacts  on  the  frame  of  the  device  through  a 
knife  edge.  The  vertical  guide  sheaves  are  attached  to  the  yoke  of  the 
load  frame,  which  deflects  as  the  torque  load  increases.  All  guide 
sheaves  can  be  raised  or  lowered  so  that  they  are  in  a horizontal  plane 
with  the  torque  wheel  at  the  start  of  a test.  The  torque  loading  sys- 
tem is  shown  in  fig.  A3. 

5.  For  this  testing  program,  the  very  rapid  rates  of  displacement 
required  that  test  data  be  recorded  electronically.  To  accomplish  this, 
a linear  variable  differential  transformer  (LVDT)  was  mounted  on  the 
load  frame  to  measure  its  deflection.  This  modification  can  be  seen  in 
fig.  Al*.  Fr.^m  the  deflection  of  the  load  frame,  the  shear  stress  in  the 
soil  can  be  computed. 


Normal  Loading  System 

6.  Details  of  the  unmodified  normal  loading  system  are  shown  in 
fig.  A3.  The  unmodified  system  consisted  of  a hanger  and  counterbal- 
anced double  lever  to  which  dead  weights  were  added.  The  load  was  trans- 
mitted to  the  torque  wheel  and  annular  piston  by  a loading  rod.  A turn- 
buckle  on  this  rod  permitted  adjustment  of  the  levers  to  a horizontal 
position. 
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T.  In  this  testing  program,  soil  specimens  were  sheared  under 
constant  volume;  l.e.,  during  shear,  the  normal  stress  was  continuously 
varied  to  counteract  any  tendency  of  the  soil  specimen  to  expand  or  con- 
tract. Operating  In  this  manner  ensured  that  no  drainage  would  occur, 
and  hence  a measure  of  undralned  shearing  resistance  was  obtained.  The 
mechanisms  of  the  constant-volume  test  In  a direct  shear  box  are  dis- 
cussed In  detail  In  reference  23.  The  basic  principles  of  the  constant- 
volume  direct  shear  test  also  apply  to  the  constant-volume  annular  shear 
test. 

8.  In  order  to  record  and  control  the  normal  load,  the  lever 
system  was  replaced  by  a 1000-lb  electronic  load  cell,  which  was  held  at 
the  lower  end  by  a rod  and  bearing  that  was  fastened  to  the  frame  of  the 
device  and  which  at  the  upper  end  was  connected  to  the  loading  rod  by  a 
turnbuckle.  The  rod  and  bearing  provided  alignment  for  the  load  cell, 
and  the  turnbuckle  provided  any  necessary  load  adjustment.  The  modified 
normal  loading  system  Is  shown  In  figs.  A6  and  AT. 

Volume  Control 


9.  The  volume  of  the  annular  shear  test  specimen  was  controlled 
by  controlling  the  thickness  of  the  specimen.  A 0.001-ln.  Ames  dial  In- 
dicator was  used  to  control  the  thickness  during  shear  within  a range  of 
+0.0001  In.  The  vertical  thickness  dial  Is  held  rigidly  In  the  dial 
frame,  which  clraps  to  the  frame  of  the  device.  The  dial  stem  rests  on 
the  top  cap,  which  Is  held  In  the  torque  plate  that  Is  rigidly  clamped 
to  the  annular  piston.  The  annular  piston  rests  on  the  test  specimen. 
With  this  assembly,  any  tendency  for  vertical  expansion  or  contraction 
of  the  soil  can  be  detected  with  the  dial  Indicator,  and  the  normal 
stress  can  be  adjusted  to  counteract  this  tendency  for  movement. 

10.  The  signals  from  the  electronic  load  cell  and  the  LVDT  were 
amplified  and  recorded  on  a strip  chart  recorder.  The  composite  testing 
apparatus  Is  shown  In  fig.  A8. 
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Fig.  A8.  Modified  «nn\ilar  shear  apparatus 


APPENDIX  B:  PREPARATION  OF  RUBBER  ANNULUS 


1.  In  order  to  obtain  a realistic  measure  of  the  resistance  at 
the  Interface  of  soil  and  rubber  for  the  wheel-soil  interaction  study,  a 
rubber  annulu'  was  fabricated  in  order  that  annular  torsion  tests  could 
be  performed. 

2.  The  rubber  used  to  mold  the  annulus  was  of  the  type  used  com- 
mercially for  retreading  truck  tires.  The  annulus  used  in  this  testing 
program  was  prepared  using  a procedure  that  simulated  the  commercial  re- 
treading operation  as  closely  as  possible. 

3.  A mold  that  would  allow  the  forming  of  the  annulus  in  its 
final  dimensions  was  constructed;  it  can  be  seen  schematically  in 

fig.  Bl.  A photograph  of  the  mold  is  presented  in  fig.  B2.  A quantity 
of  uncured  rubber  available  commercially  from  the  General  Tire  and  Rub- 
ber Company  was  obtained  and  cut  into  an  annular  geometry,  1/2  in.  thick, 
which  would  fit  into  the  mold.  The  rubber  was  bonded  to  the  cleated 
stainless  steel  plate  (fig.  B3).  The  rigid  plate  ensured  that  the  rub- 
ber annulus  would  keep  its  geometric  shape  and  provided  a means  of  me- 
chanical linkage  to  the  shear  apparatus. 

t«.  The  steel  plate  was  put  into  the  bottom  of  the  mold,  the  annu- 
lus of  uncured  rubber  was  put  in  on  top  of  the  plate,  and  a smooth  steel 
ring  was  put  on  top  of  the  uncured  rubber.  The  mold  was  then  completely 
assembled  and  a torque  applied  to  the  top  nut  of  the  assembly  to  give  a 
pressure  on  the  rubber  of  l60  psi.  The  assembly  was  then  placed  in  an 
oven  at  1U9  C (or  about  300  F)  and  allowed  to  cure  at  this  temperature 
for  a period  of  U3  min.  The  assembly  was  then  allowed  to  cool,  and  the 
cured  rubber  annulus  and  its  steel  linkage  plate  forced  from  the  mold. 

The  final  product  is  shown  in  fig.  Bi». 
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Finished  nihber  annulus  shown 


